ADA105624 


AFWAL-TR-80-31 36 


VALUES  AND  VECTORS  FROM  SINUSOIDAL 
EXCITATIONS  AT  NEAR  RESONANCE  FREQUENCIES 


C.  L.  Keller 

Applied  Mathematics  Group 
Analysis  and  Optimization  Branch 


June  1981 


Interim  Report  for  Period  August  1978  -  June  1980 


Approved  for  public  release;  distribution  unlimited. 


FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND  ,, 

WRIGHT- PATTERSON  AIR  FORCE  BASE,  OHlQ4M33 


14 


NOTICE 


when  Government  drawings,  specifications ,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
anu  way  supplied  the  said  drawings,  specifications ,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS ,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  techrdcal  report  has  been  reviewed  and  is  approved  for  publication. 


<J2  A 


CHARLES  L.  KELLER 
Mathematician 
Applied  Mathematics  Group 
Analysis  &  Optimization  Branch 


FOR  THE  COVMANDER 


RALPH  L.  KUSTER,  Jr.,  Colonel,  USAF 
Chief,  Structures  &  Dynamics  Division 


/7? 


FREDERICK  A.  PICCHIONI,  Lt  Col,  USAF 
Chief,  Analysis  &  Optimization  Branch 
Structures  &  Dynamics  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  you r  organization  please  notift/AFWAL/F_IBRD» 
W-PAFB,  OH  45431  to  help  us  maintain  a  current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


<it'  IJPI  r  /  •  1  ASSlFT  A  TIOll  O F  THIS  PAGE  flWion  D»/»  Entered) 

1 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

| 

2 

I  REPORT  NUMBER  j.  GOVT  ACCESSION  no. 

SAfTJftnfR-8tf-3136  j  CDAIDSGxZ.1) 

3.  RECIPIENT’S  CATALOG  NUMBER 

/ 

4A  TITLE  (end  Subtitle)  S~ 

]  DETERMINATION  OF  COMPLEX  .CHARACTERISTIC:  (7, 

-'f  VALUES  AND  ACTORS  FROM  SINUSOIDAL  .-*“77* 

ii  TT°g  ~*  RyPOWT-B- PERIOD  COVERED 

Interim  ^pC’t*  ' - - 

WMiv-  78  —  June  80,,] 

j  EXCITATIONS  AT  NEAR  RESONANCE  (ft 

J  Irequencies#  ,  “* 

7.  AuTMO«r») 

C.  L.J  Keller/ 

9.  contract  or  grant  number^; 

In-House  Research 

9.  performing  organization  name  and  address 

Flight  Dynamics  Laboratory  (AFWAL/FIBRD) 

Air  Force  Wright  Aeronautical  Laboratories  (AFSC) 
Wriqht-Patterson  Air  Force  Base,  Ohio  45433 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 

AREA  ft  WORK  UNIT  NUMBERS  _ _ _ 

DoD  Element  61102F 
r?7OTH!Q2 - - - 

U.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Flight  Dynamics  Laboratory  (AFWAL/FIBR)  Q j 

Air  Force  Wright  Aeronautical  Laboratories  (AFSut'' 
Wriqht-Patterson  Air  Force  Base,  Ohio  45433 

9SX&  rr 

13.  NUMBER  OF  PAGES  /J,  1  J'  * 

133  Her  A  -  .  -  - 

1«.  MONITORING  AGENCY  NAME  ft  AOORESSf/l  different  Irom  Controlling  Ollier,) 

15.  SECURITY  CLASS,  (olthle  report) - 

UNCLASSIFIED 

is*  declassification/oowngraoing 
SCHEDULE 

16  DISTRIBUTION  STATEMENT  (ot  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ot  the  abetted  entered  In  Block  20,  It  dlllerent  Iron,  Report) 

19.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  revert*  aid*  1/  neceaaery  and  Identity  by  block  number) 

Modal  Parameters 

Linear  Systems  Identification 

Mass,  Stiffness  and  Damping  Matrices 

Ground  Vibration  Testing 

Dynamic  Structural  Testing 

20.  ABSTRACT  (Continue  on  reverie  elde  It  neceeeery  end  Identity  by  block  number) 

"~Tn  this  report  an  iterative  procedure  for  determining  complex  characteristic 
values  and  vectors  from  measured  values  of  the  steady  state  response  to  sinus¬ 
oidal  excitations  is  described.  An  extension  of  the  procedure  which  enables  it 
to  cope  with  the  case  of  two  close  characteristic  values  is  presented  also. 
Results  of  numerical  experiments  performed  to  investigate  these  procedures  are 
given.  This  report  contains  also  extensive  background  material  and  an  exami¬ 
nation  of  the  methods  of  Ibrahim,  of  Wittmeyer  and  of  Link  and  Vollan. 

....  4 - 

DD  .  1473  EOITION  OF  1  NOV  «S  IS  OBSOLETE  '  ^  \ 

DO  ,  :r,3  1473  EOITION  OF  I  NOV  •*  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWi»n  Dele  Entered) 


AFWAL-TR-80-31 36 


FOREWORD 

This  report  describes  work  performed  in  the  Applied  Mathematics  Group 
of  the  Analysis  and  Optimization  Branch,  Structures  and  Dynamics 
Division  of  the  Flight  Dynamics  Laboratory  (AFWAL/FIBRD)  under  Project 
2304N1  Computational  Aspects  of  Fluid  and  Structural  Mechanics,  Work  Unit 
2304N102.  This  is  an  interim  report  on  work  carried  out  between 
August  1978  and  June  1980.  The  author,  C.  L.  Keller,  submitted  the 
report  in  October  1980. 

The  author  wishes  to  thank  Dr.  K.  G.  Guderley  for  many  stimulating 
discussions  and  suggestions.  He  thanks  also  Ms.  Mary  Li  pi k  for  her  . 
typing  and  assistance  in  preparing  this  report. 


\ 


ii  i 


fflSCBttKl  PAOS  BUNK-MOT  71JLMKD 


1 


AFWAL-TR-80-31 36 


TABLE  OF  CONTENTS 

SECTION  PAGE 

I  INTRODUCTION  1 

II  DETERMINATION  OF  THE  CHARACTERISTIC  VALUES  AND  VECTORS  7 

III  TWO  CLOSE  CHARACTERISTIC  VALUES  17 

IV  THE  FREQUENCY  RESPONSE  FUNCTION  BY  NEWTON’S  METHOD  21 

V  RESULTS  AND  CONCLUSIONS  25 

APPENDIX  A  LINEAR  SYSTEMS  OF  ORDINARY  DIFFERENTIAL 

EQUATIONS  29 

APPENDIX  B  INTERPOLATION  WITH  EXPONENTIAL  FUNCTIONS  42 

APPENDIX  C  SOME  MATRIX  EIGENVALUE-EIGENVECTOR 

CONSIDERATIONS  53 

APPENDIX  D  DETERMINATION  OF  MASS,  DAMPING  AND 

STIFFNESS  MATRICES  59 

REFERENCES  123 


v 


mCUNB  PAOS  BUML-NOT  FILMED 


AFWAL-TR-80-31 36 


LIST  OF  ILLUSTRATIONS 


FIGURE 

PAGE 

1 

An  Example  of  an  Element  of  the  Frequency  Response 

Function  Matrix 

68 

2 

Graph  of  the  Magnitude  of  w(w) 

69 

3 

Graph  of  the  Argument  of  w(u>) 

70 

4 

Graph  in  the  Complex  Plane  of  t(u)  +  inU)  =  w(w)  -  b4/(iw-X4) 

70 

5 

Graph  of  |  w(to)  | 

and  of  |w(u>)  -  b^/ ( ito-X^)  | 

71 

6a 

Graphs  of  w(w)  for  Four  Values  of  b 

4 

72 

6b 

Graphs  of  |w(oj) 

|  for  the  Functions  w(u)  of  Figure  6a 

73 

7a 

Graph  of  w(aO , 

b4  =  0.866025404  +  0. 5i 

77 

7b 

Graph  of  | w (to)  | 

,  b4  =  0.866025404  +  0.5i 

78 

8a 

Graph  of  w(w). 

b4  =  0.866025404  -  0.5i 

82 

8b 

Graph  of  |w(co)  | 

,  b4  =  0.866025404  -  0.5i 

83 

9a 

Graph  of  w (w). 

b4  =  -0.866025404  -  0.5i 

87 

9b 

Graph  of  j w{w) | 

,  b4  =  -0.866025404  -  0.5i 

88 

1  Oa 

Graph  of  w(io) , 

b4  =  -0.866025404  +  0.5i 

92 

10b 

Graph  of  |w(u>)  | 

,  b4  =  -0.866025404  +  0.5i 

93 

11a 

w{u>)  for  Four  Values  of  Xg 

94 

11b 

Graphs  of  (w(oj) 

|  for  the  Functions  w(w)  of  Figure  11a 

95 

1  ?a 

Graph  of  w(w) , 

Xg  =  -0.1  +  4.5i 

99 

12b 

Graph  of  |w(w) | 

,  Ag  =  -0.1  +  4.5i 

100 

13a 

Graph  of  w(w). 

Ag  =  -0.1  +  4. 3i 

104 

13b 

Graph  of  jw(u>)  | 

,  Ag  =  -0.1  +  4 . 3i 

105 

14a 

Graph  of  w(u) , 

Ag  =  -0.1  +  4.2i 

109 

14b 

Graph  of  |w(w)| 

,  Ag  =  -0.1  +  4.2i 

110 

vi  i 


HOCUS*?  p AOS  BUMC-N0T  71 


AFWAL-TR-80-3136 


LIST  OF  ILLUSTRATIONS  (Concluded) 

FIGURE 

PAGE 

15a 

Graph  of  w(co),  Xg  =  -0.1  +  4.11 

114 

15b 

Graph  of  |w(oj)|,  Xg  =  -0.1  +  4.11 

115 

16a 

Graph  of  w(oj),  Xg  =  -0.1  +  4.05i 

119 

16b 

Graph  of  |w(u)|,  Xg  =  -0.1  +  4.05i 

120 

17 

Comparison  of  Values  of  w(<o)  for  w(ui)  Computed  Using 
the  Exact  Characteristic  Values  and  Coefficient  and  w(w) 
Computed  Using  the  Characteristic  Values  and  Coefficients 
from  TABLE  9C 

122 

! 


viii 


V 


AFWAL-TR-80-31 36 


LIST  OF  TABLES 

TABLE  PAGE 

1A  Defining  Parameters,  Function  Values  and  Truncated 

Values  t>4  =  0.86602E404  +  0.51  74 

IB  Computed  Defining  Parameters  from  "Exact"  Function 

Values  75 

1C  Computed  Defining  Parameters  from  Truncated  Function 

Values  76 

2A  Defining  Parameters,  Function  Values  and  Truncated 

Values  =  0.866025404  +  0.5i  79 

28  Computed  Defining  Parameters  from  "Exact"  Function 

Values  80 

2C  Computed  Defining  Parameters  from  Truncated  Function 

Values  81 

^  3A  Defining  Parameters,  Function  Values  and  Truncated 

Values  t>4  =  -0.866025404  -  0.5i  84 

3B  Computed  Defining  Parameters  from  "Exact"  Function 

*  '  Values  85 

^  3C  Computed  Defining  Parameters  from  Truncated  Function 

Values  86 

4A  Defining  Parameters,  Function  Values  and  Truncated 

Values  =  -0.866025404  +  0.51  89 

4B  Computed  Defining  Parameters  from  "Exact"  Function 

Values  90 

4C  Computed  Defining  Parameters  from  Truncated  Function 

Values  91 

5A  Defining  Parameters,  Function  Values  and  Truncated 

Values  Ag  ~  -0.1  +  4.5i  96 

5B  Computed  Defining  Parameters  from  "Exact"  Function 

Values  97 

5C  Computed  Defining  Parameters  from  Truncated  Function 

Values  98 

6A  Defining  Parameters,  Function  Values  and  Truncated 
>  Values  Ag  =  -0.1  +  4.3i  101 

68  Computed  Defining  Parameters  from  "Exact"  Function 

Values  102 


V 


LIST  OF  TABLES  (Concluded) 


TABLE 

6C  Computed  Defining  Parameters  from  Truncated  Function 
Val ues 

7A  Defining  Parameters,  Function  Values  and  Truncated 
Values  Ag  =  -0.1  +  4.2i 

7B  Computed  Defining  Parameters  from  "Exact"  Function 
Values 

7C  Computed  Defining  Parameters  from  Truncated  Function 
Values 

8A  Defining  Parameters,  Function  Values  and  Truncated 
Values  Ag  =  -0.1  +  4.1  i 

8B  Computed  Defining  Parameters  from  "Exact"  Function 
Val ues 

8C  Computed  Defining  Parameters  from  Truncated  Function 
Val ues 

9A  Defining  Parameters,  Function  Values  and  Truncated 
Values  Ag  =  -0.1  +  4.05i 

9B  Computed  Defining  Parameters  from  "Exact"  Function 
Values 

9C  Computed  Defining  Parameters  from  Truncated  Function 
Values 

9D  Function  Values  Determined  Using  the  Computed 
Characteristic  Values  and  Coefficients  From 
TABLE  9C 


AFWAL - TR- 80-31 36 


aagarar 


SECTION  i 
INTRODUCTION 

In  the  last  10  years  a  number  of  procedures  have  been  proposed  for 
determining  vibration  parameters  of  a  structure  from  its  measured 
response  to  known  excitations.  The  reason  for  this  activity  was  due  to 
shortcomings  and  dissatisfaction  with  the  classical  phase  resonance  or 
tune  and  dwell  procedure. 

In  the  past,  the  phase  resonance  method  was  characterized  as  requiring 
a  rather  lengthy  set  up  time.  Practical  limits  on  the  number  and 
positioning  of  the  shakers  result  in  limitations  on  the  accuracy  of  the 
method.  The  desired  information  is  obtained  directly  from  the  test. 
However,  obtaining  satisfactory  test  results  was  a  slow  process  requiring 
highly  skilled  operators.  Consequently,  the  test  structure  was  not 
available  for  other  programed  activities  for  a  lengthy  period,  usually, 
at  a  crucial  time. 

Accordingly,  other  methods  were  sought  for  which  the  test  set  up  would 
be  simpler,  require  fewer  excitation  points,  and  be  less  dependent  upon 
operator  skill.  The  desired  information  need  not  be  given  directly  by 
the  test  but  would  be  obtained  at  a  later  date  from  the  test  data  after 
the  structure  was  released. 

For  those  methods  for  which  test  results  are  not  known  until  some 
time  after  the  structure  is  released  there  is  a  danger  that,  due  to  an 
undetected  error  in  the  test  set  up,  the  test  data  is  worthless. 
Frequently,  in  such  cases  it  is  impossible  to  reschedule  the  test  until 
after  such  a  time  that  the  main  reason  for  the  test  is  past. 

In  Section  II  of  this  report  a  procedure  is  described  for  computing 
the  complex  characteristic  values  and  vectors  from  the  steady  state 
response  to  sinusoidal  excitations  at  a  single  point.  In  Section  III 
we  describe  a  modification  of  the  procedure  so  that  it  can  handle  the 
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case  of  two  close  characteristic  values.  In  a  previous  report, 
AFFDL-TR-78-59  (Reference  1),  we  have  shown  how  to  modify  the  procedure 
to  handle  characteristic  values  of  multiplicity  greater  than  1. 

We  believe  the  procedure  described  in  Sections  II  and  III  has  most  of 
the  advantages  of  the  phase  resonance  method  with  most  of  the  shortcomings 
eliminated.  Since  it  depends  upon  excitations  at  a  single  point  the  test 
set  up  is  simpler.  Also,  since  data  only  in  the  vicinity  of  resonance  is 
used,  only  the  frequency  needs  to  be  adjusted.  The  frequency  adjustments 
do  not  have  to  be  too  precise,  consequently,  the  degree  of  operator  skill 
is  of  less  importance. 

The  computations  are  simple.  We  believe  the  procedure  can  be 
automated  so  that  the  final  results  are  obtained  practically  on  line. 
Accordingly,  errors  in  test  set  up,  if  any,  should  be  detected  early  in 
the  allotted  time  for  the  tests,  corrected  and  the  test  rerun.  Moreover, 
the  nonlinear  damping  and  stiffness  properties  can  be  calculated  as 
functions  of  the  excitation  amplitudes. 

The  procedure  described  in  Sections  II  and  III  is  based  upon  exci¬ 
tations  at  a  single  point.  However,  it  should  be  clear  from  the 
discussion  of  these  sections  that  the  procedure  is  not  limited  to  single 

point  excitation.  For  multipoint  excitation,  the  quantity  v^r,  Equation  2, 

K  T 

no  longer  isolates  a  single  component  of  the  vector  v^  but  instead,  v^r 
is  some  linear  combination  of  the  components  of  v^.  Thus  for  multipoint 
excitation  one  needs  to  perform  sufficiently  many  additional  independent 
experiments  to  determine  at  least  one  (and  consequently  all)  of  the 
components  of  v^  appearing  in  the  linear  combination  v^r.  It  is  clear 
also  that  this  is  accomplished  by  varying  components  of  the  vector  r,  v^r 
that  is,  by  varying  the  amplitude  at  the  set  of  multiple  excitation  points 
in  such  a  way  that  the  total  energy  input  is  essentially  constant. 

Using  the  procedure  of  Sections  II  and  III  modified  for  multiple 
point  excitation  does  not  affect  the  test  set  up  appreciably.  Rather,  it 
amounts  to  some  additional  experiments  and  computations  which  are  essen¬ 
tially  of  the  same  character  as  for  single  point  excitation. 
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In  Section  IV  we  describe  a  Newton  procedure  for  determining  the 
characteristic  values  and  make  some  additional  remarks  concerning  the 
iteration  process  of  Section  II.  In  Section  V  we  describe  the  numerical 
experiments  which  we  perform  to  test  the  procedure  of  Sections  II  and  III 
and  give  the  results  of  these  numerical  experiments. 

In  Appendices  A-D  we  give  mathematical  background  for  the  method  of 
Section  II  and  for  some  of  the  other  methods  examined.  Appendix  A  is 
mainly  results,  from  the  theory  of  linear  systems  of  differential 
equations  with  constant  coefficients,  which  are  helpful  for  the  identi¬ 
fication  problem.  Appendix  B  develops,  in  considerable  detail,  Prony's 
method  and  contrasts  Prony's  method  with  the  method  of  Ibrahim 
(References  2-5)  for  determining  the  parameters  in  a  sum  of  exponential 
functions.  Appendix  C  considers  some  aspects  of  the  matrix  eigenvalue- 
eigenvector  problem  and  the  use  thereof  in  Wittmeyer's  method 
(References  6  and  7).  In  Appendix  D  we  discuss  the  problem  of  solving 
systems  of  linear  equations  which  may  be  over,  even,  or  under  determined 
and  the  method  Link  and  Vollan  (Reference  8). 

The  methods  for  determining  vibration  parameters  can  be  classified 
in  various  ways  depending  upon  which  features  are  emphasized.  One 
rather  evident  category  consists  of  those  methods  which  determine  mode 
shapes  and  frequencies.  That  is,  numbers  A  and  vectors  u  which  satisfy 
the  conditions,  either 

[A2M  +  AC  +  K]u  =  0 


or 

[A2M  +  K]u  =  0 

A  second  category  consists  of  those  methods  which  determine  the  mass, 
damping  and  stiffness  matrices  M,  C  and  K  directly.  The  second  category 
appears  to  consist  of  variants  of  the  Link  and  Vollan  method  and  the 
"ill-conditioned"  method  described  in  Appendix  D.  Thus  most  methods 
belong  to  the  first  category. 
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Many  methods  use  the  frequency  response  function  (Equation  A38)  in 
some  way  or  other.  The  frequency  response  function  can  be  written  in 
different  ways  depending  on  how  the  characteristic  vectors  are  normalized 
and  which  vibration  parameters  are  displayed. 

The  method  discussed  in  References  9-13  computes  (a  column  of)  the 
frequency  response  function  matrix  from  the  response  to  excitation  of 
impulse  or  random  type.  The  procedure  uses  fast  Fourier  analysis 
methods  (Equation  A28)  and  is  highly  automated.  In  this  procedure  the 
frequency  response  function  is  obtained  numerically,  that  is,  as  a  set 
of  number  pairs.  The  second  entry  of  this  number  pair  is  the  "value" 
of  the  frequency  response  function  and  the  first  entry  is  the  frequency 
to  to  which  the  second  entry  corresponds.  Lastly,  the  parameters  in  the 
frequency  response  function  are  determined,  usually,  according  to  the 
least  squares  criteria.  Thus  advances  in  electronics,  in  sensing  and 
recording  equipment  and  minicomputers,  have  made  a  complicated  procedure 
economically  practical. 

Hence  the  fast  Fourier  analysis  procedure  discussed  in  References  9-13 
uses  broad  band  excitation  to  produce  a  transient  response.  From  the 
Fourier  analysis  of  this  experimentally  obtained  data  one  obtains  the 
frequency  response  function  (numerically).  Finally,  from  the  "fit" 
of  the  expression  for  the  frequency  response  function  to  the  experimentally 
determined  values  of  the  frequency  response  function  one  obtains  the 
complex  mode  shapes  and  frequencies. 

The  method  of  Section  II  is  based  upon  the  frequency  response 
function  matrix  and  its  properties  also.  The  complex  amplitude  y  of 
the  steady  state  response  y  exp  (iwt)  to  the  harmonic  excitation  r  exp 
(iwt)  is,  essentially,  the  value  of  the  frequency  response  function 
corresponding  to  the  frequency  u>.  From  knowledge  of  the  nature  of  the 
frequency  response  function  we  choose  a  set  of  frequencies  u>  and  obtain 
the  values  of  the  frequency  response  function  corresponding  to  the  set 
of  frequencies  w.  These  values  of  the  frequency  response  function  are 
obtained  experimentally  as  the  complex  amplitude  y  of  the  steady  state 
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response  to  the  harmonic  excitation  r  exp  (iwt).  The  set  of  values  of 
the  frequency  response  function,  so  obtained,  are  used  to  set  up  a 
system  of  nonlinear  equations  for  the  parameters  (coefficients)  of  the 
frequency  response  function.  The  complex  mode  shapes  and  frequencies 
are  obtained  as  the  solution  to  this  nonlinear  system  of  equations. 

The  method  of  Wittmeyer  (References  6  and  7)  is  yet  another  method 
based  upon  the  frequency  response  function  and  its  properties.  This 
procedure  is  a  method  for  solving  a  matrix  eigenvector-eigenvalue  problem 
appropriately  modified  to  be  applicable  to  the  ground  vibration  test 
problem.  Matrix  eigenvector-eigenvalue  methods  are  iterative  procedures, 
generally.  The  experimental  data  is  obtained  from  multipoint  sinusoidal 
excitation.  One  alternates  between  experiment  and  computation.  The 
current  Approximation  to  the  eigenvector  (being  determined)  is  obtained 
from  experiment.  If  this  approximation  is  not  satisfactory  an  improved 
excitation  vector  is  computed  from  the  experimental  data  and  the 
experiment  repeated. 

The  method  of  Link  and  Vollan  (Reference  8)  determines  the  damping 
and  stiffness  matrices  directly  from  steady  state  response  to  sinusoidal 
excitation  data.  (The  mass  matrix  M  usually  is  assumed  known.)  This 
procedure  uses  the  same  kind  of  data  as  the  method  of  this  report.  The 
damping  and  stiffness  matrices  C  and  K  are  obtained  as  the  solution  to 
a  system  of  linear  equations. 

Using  the  pseudo  or  generalized  inverse,  if  necessary,  the  linear 
system  of  equations  is  solvable  for  the  matrices  C  and  K  regardless  of 
whether  the  system  of  equations  for  C  and  K  is  over,  even  or  under 
determined.  Of  course,  if  the  system  is  under  determined  then  the 
matrices  C  and  K  are  even  less  well  known. 

An  adequate  number  of  linearly  independent  conditions  are  needed 
for  a  strict  solution  for  the  matrices  C  and  K.  Equation  A29  shows  the 
relation  of  the  matrices  M,  C  and  K  to  the  frequency  response  function. 
Equation  A29  also  shows  that  linearly  independent  conditions  are  obtained 
at,  or  near  to  resonance  frequencies. 
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Ibrahim's  method  (References  2-5)  determines  the  complex  eigenvectors 
and  frequencies  from  the  free  response  of  the  system.  Two  matrices  are 
generated  from  the  data  obtained  from  the  free  response  sampled  at  equal 
time  increments.  The  complex  mode  shapes  and  frequencies  are  obtained 
from  the  eigenvalue  and  eigenvector  solution  to  the  generalized 
eigenvalue-eigenvector  problem  associated  with  these  two  matrices. 

The  method  produces  an  equivalent  model  for  the  structure.  That  is, 
it  produces  a  model  which  has  the  same  eigenvectors  as  the  structure. 
However,  since  no  particular  normalization  is  required  or  enforced,  the 
mass,  stiffness  and  damping  matrices  cannot  be  determined  by  the  method; 
at  least,  not  without  some  additional  experiments  and  computation. 

The  phase  resonance  method  has  also  made  use  of  the  development  in 
computers  and  visual  display  techniques  to  bring  about  improvements  in 
the  method.  The  indicator  function,  basically,  is  the  sun  of  the 
magnitudes  of  all  the  in  phase  amplitudes  at  all  the  measurement  points. 
With  the  on  line  computation  and  display  of  the  indicator  function  an 
operator  can  see  immediately  the  effect  of  a  change  in  amplitude  at  an 
excitation  point  or  the  effect  of  a  change  in  the  excitation  frequency. 
The  phase  resonance  method  is  still  widely  used  in  ground  vibration 
tests. 
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SECTION  II 

DETERMINATION  OF  THE  CHARACTERISTIC  VALUES  AND  VECTORS 

In  this  section  a  method  is  described  for  determining  the  parameters 
in  an  expression  for  the  frequency  response  function  matrix  from  the 
steady  state  response  to  sinusoidal  excitations.  The  quantities 
determined  are  the  characteristic  values  A^  k=l ,  m,  the  corresponding 

characteristic  vectors  u^  and  the  corresponding  characteristic  vectors 
v^  of  the  transposed  system.  (The  notation,  conventions  and  results  of 
Appendix  A  are  assumed  as  being  familiar  throughout  this  report.) 

The  method  described  here  is  essentially  the  same  as  in  Reference  1. 
However,  in  Reference  1  vectors  were  treated  as  single  entities.  That 
is,  very  little  attention  was  paid  to  the  individual  components  except 
for  two  instances  when  the  fact  was  used  that  one  could  assume  a  special 
value  or  special  values  for  a  component  or  several  components.  Here,  in 
order  to  display  certain  features  with  clarity,  the  process  is  described 
in  terms  of  a  scalar  component. 

The  frequency  response  function  parameters  will  be  determined  by  an 
iterative  process.  We  do  not  have  a  formal  proof  for  the  convergence  of 
this  process.  However,  we  are  able  to  make  some  observations  which  lead 
one  to  believe  that  the  process  is  convergent.  Additionally,  our 
numerical  experience  with  the  iteration  process,  although  limited,  has 
not  indicated  any  difficulties,  thus  providing  evidence  that  the  process 
is  convergent  under  conditions  sufficiently  general  for  our  purposes. 

Usually,  in  determining  the  characteristic  values  and  vectors,  one 
assumes  that  the  characteristic  values  are  well-separated.  In  this 
section  we  make  this  assumption  also.  In  Reference  1  we  showed  that  the 
method  described  here  could  be  modified  to  handle  characteristic  values 
of  multiplicity  greater  than  1.  In  Section  III  of  this  report,  we  will 
show  how  to  extend  the  procedure  of  this  section  so  that  it  can  handle 
the  case  of  two  close  characteristic  values  also. 
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The  procedure  described  in  References  9-13  analyzes  the  transient 
structural  response  to  broadband  excitation  and  uses  the  relation  between 
the  frequency  response  (or  transfer)  function  matrix  and  the  transforms 
of  the  excitation  and  response  as  expressed  by  Equation  A28.  The 
procedure  to  be  described  in  this  section  uses  instead  the  relation 
between  the  frequency  response  function  matrix,  the  steady  state  response 
to  an  harmonic  excitation  and  the  harmonic  excitation. 

It  is  assumed  that  after  discretization  the  equations  of  motion  of 
the  structure  are  adequately  represented  by  a  system  of  equations  of 
the  form 

Mx  +  Cx  +  kx  =  f  (1 ) 

Here  M,  C  and  K  denote  the  mass,  damping  and  stiffness  matrices 
respectively.  M,  C  and  K  are  taken  as  square  matrices  of  order  m.  In 
addition,  these  matrices  are  usually  assumed  to  be  symmetric.  The 
components  of  the  vector  x  usually  are  displacements  at  stations  on  the 
structure . 


The  steady  state  response  of  a  system  of  equations  of  the  form  of 
Equation  1  to  an  harmonic  excitation  f(t)=r  exp(icot)  is  x(t)=y  exp(iu)t). 
Here  r  and  y  denote  constant  vectors.  The  frequency  response  function 
matrix  and  the  vectors  r  and  y  satisfy  the  condition 


;  Vkr 

y  =  l  - r- 

k=l  w  ‘  Xk 


(2) 


(Equation  A38).  The  vector  y  is  determined  experimentally  in  some  way  from 
the  steady  state  response  to  the  excitation  r  sin  cat.  One  possible  way 
is  from  the  system  of  equations  (Equation  A37).  The  vector  r,  of  course, 
is  chosen,  usually  with  only  one  component  differing  from  zero.  This 
is  the  component  corresponding  to  the  station  on  the  structure  at  which 
the  excitation  is  applied. 


The  frequency  response  function  matrix  may  be  written  in  various  ways. 
The  particular  form  is  a  consequence  of  the  normality  conditions  imposed 
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on  and  satisfied  by  the  characteristic  vectors  and  vk,  k=l ,  . ..,  m. 
For  characteristic  vectors  which  satisfy  the  normality  condition 
expressed  by  Equation  A13,  the  frequency  response  function  matrix  takes 
the  form  given  in  Equation  2.  It  should  be  clear  that  one  does  not  need 
to  know  the  vectors  u^  and  v^  normalized  so  that  they  satisfy  the 
condition  (Equation  A13).  Rather  u^  and  v^  are  determined  so  that 
Equation  2  is  satisfied  and  when  so  determined,  then  the  normality 
conditions  (Equation  A13)  are  satisfied  also. 


This  enforcement  of  the  normality  conditions  enables  one  to  determine 
the  matrices  M,  C,  and  K  after  the  quantities  X^,  u^  and  v^  are 
determined.  Thus  after  X^,  u^  and  v^  are  determined  so  that  Equation  2 
is  satisfied  for  the  steady  state  response  to  any  harmonic  excitation, 
the  relations  expressed  by  Equations  A22,  A24  and  A25  hold.  It  should 
not  be  inferred  that  these  relations  will  be  used  directly  in  computing 
M,  C  and  K.  Rather  some  relations  derived  from  Equations  A22,  A24  and 
A25  or  the  procedure  given  in  References  12  and  1  may  be  more  efficient 
for  computing  these  matrices. 


From  Equation  2  it  is  clear  that  for  a  vector  r  with  only  the  pth 

component  different  from  zero,  the  vector  y  is  this  pth  component  times 

the  pth  column  of  the  frequency  response  function  matrix.  It  follows 

then  that  the  jth  component  of  the  vector  y  is  a  complex  valued  function 

w.(u>)  of  the  form 
J 


Wj  (w) 


m  m 

E  b  . /(icu-X  )  +  E 

k=l  K  k=l 


V(1u‘V 


(3) 


If  the  b^  and  X^  are  not  known,  then  the  expression  for  the  function 

Wj(w)  contains  n=2m  unknowns.  To  determine  the  b^  and  X^,  n  independent 

conditions  are  needed.  These  n  conditions  are  provided  by  the  values  of 

Wj(cu)  corresponding  to  n  distinct  values  of  w.  We  will  show  that  if  we 

know  w.(w)  for  n  judiciously  chosen  values  w  =  to  ,  p=l ,  ...,  n  the 
J  P 

resulting  system  of  nonlinear  equations  is  readily  solved  for  the  X^ 

and  b^  by  an  iteration  process. 
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i 


In  Figure  1  we  have  plotted  in  a  complex  plane  values  of  a  function 
w(w)  of  the  form  expressed  by  Equation  3  for  values  of  to  ranging  from  0.5 
to  12  at  increments  of  0.05.  Figures  2  and  3  are  plots  of  values  of  the 
magnitude  and  argument  respectively  of  the  same  function  for  the  same 
range  of  values  of  to.  These  figures  display  the  features  which  charac¬ 
terize  functions  w.(to)  of  the  form  expressed  by  Equation  3.  The 
3 

appearance  of  Figure  3  is  determined  by  the  range  of  values  selected 
for  the  argument  of  a  complex  number. 

From  Figure  2  it  is  observed  that  the  local  maximum  of  |w(co)|  occurs 

at  some  value  of  to  close  to  Im[ A^]  for  k=l ,  ....  m.  From  Equation  A36  it 

follows  that  |w.(co)|  is  the  amplitude  of  the  steady  state  response  to  r 
3 

sin  tot.  Also,  in  Figure  1,  if  we  would  label  the  plotted  points  with 
the  value  of  to  to  which  they  correspond  -  we  would  observe  that  w(to) 
changes  most  rapidly  when  to  is  close  to  Im[A^].  That  is,  in  Figure  1 
the  values  of  w(to)  for  which  to  is  close  to  Im[X^]  are  those  points  which 
are  connected  by  rather  large  straight  line  segments. 


From  the  above  observations  it  is  clear  that  we  may  suppose  that  n 

values  of  to  and  the  corresponding  vector  y  are  known  in  pairs,  say  tOp 

and  to  ,  ,  with  corresponding  vectors  y  and  y  ,  ,  for  p=l  ,  ....  m. 

m  +  p  y  a  Jp  ■7m  +  p  k  >  » 

The  values  of  to  satisfy  the  conditions 
P 


to 


ImOp] 


<  to 


m  +  p 


(4) 


And  the  difference  to  t  -to  is  small  relative  to  the  differences 

m  +  p  p 

u>  ,  i  -  to  .  and  to  -  to  „  .  Hence  the  value  of  the  complex 

p+lm  +  p  pm  +  p-  i 

function  w^io),  which  represents  the  jth  component  of  the  vector  y(io), 
is  known  for  n  values  of  to. 


It  was  noted  above  that  the  system  of  equations  obtained  from  Equation  3 
when  w(<jo)  is  known  for  n  values  of  to  is  nonlinear.  We  will  now  describe 
an  iterative  process  for  determining  values  of  A k  and  b^k  which  satisfy 
the  system.  Set 


^  m  — 

V«)  =  kZ1  ^bjk^ia,'V  +  kjk^(iu'\))  +  &jp/(iw-Ap)  (5) 

Mp 
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Then 


bjp/(it)'V  =  Wj(a))  •  "jp<“> 


(6) 


Solving  Equation  6  for  b^  we  obtain  therefrom  the  equality 

<iv*p)[„.Up)  -  wJpUp)]  -  (iVp-»p>t«j(«Wp)  -  Vvp>!  <7> 

These  equations,  Equation  7  solved  for  Ap  and  Equation  6  solved  for  b^ 
are  the  basic  equations  of  the  iteration  process. 

We  show  next  how  to  determine  starting  values  or  initial  approximations 

for  A  and  b.  .  The  principal  observation  is  that  if  o>  and  w  ,  are 
P  JP  p  m  +  p 

close  to  Im[Ap]  then,  from  Equation  3 


and 


Wj(w  )  =  bjp/(i‘i)p-^p)  approximately 


wj(“m+p)  =  bjp/(iVp-y  approximately. 


Hence 


\(,IWV/(1VV 


which  is  readily  solved  for  Ap,  see  Equation  10  below.  Alternatively, 
from  Equation  7  we  obtain 


y"j<Vp>  -  •  'VW  ‘  "jp(“p,)]  ■ 


VWJ(lW  -  wj("p>  -  (fijp(“mtp>  -  "jp("p),]  (8) 


*  1(“«p  -  "p)[wj(‘Vp)  ■  V'Vp’1 


11 


From  the  graph  of  w^w),  Figures  4  and  5,  it  is  evident  that 

w.  (w  .  )  -  w .  (w  )  is  small,  especially  if  -  w  is  small.  Hence 

JP  m+p'  jp'  P  A  m  p  p 

neglecting  the  difference  w.  (10  )  -  w.  (cj  ),  Equation  8  can  be 

Jr  r  J  r  r 

written  as 

'  "j("p,]  "  VW  -  ViS1  ’  ,("mV“p)“jP<‘W(9' 


Also,  i f  o>  -  co  is  small  and  w.  (u>  ,  )  is  not  too  large,  we  may 

m+pp  jp  m  +  p 

neglect  the  last  term  of  Equation  9  and  we  have  as  the  equation  for  an 
initial  value  for  A 

P 

xp  *  1[“ratpwj(“m*p)  '  “p"i(“p)1/["j<‘W  '  WJ(“p)]  00) 

A 

Next,  neglecting  the  term  wjp(°J)  in  Equation  6  the  starting  value  for  b^p 
is  given  by 

bjp  ■ ‘VWV  (,1) 

Thus,  from  Equations  10  and  11  we  can  determine  initial  approximations 

A  and  b.  for  p=l ,  m  and  for  j-1 ,  ...,  m.  These  initial  values 

P  JP 

depend  on  the  function  Wj(to)  at  n  points  u>p  and  +  p  for  p=l  ,  m. 

Once  these  initial  values  are  determined,  they  can  be  used  in  the  right 
hand  side  of  the  exact  equations,  namely 

,  iW>V“mV-"jp(“W>  -  "P!wj(VjJjp(V>2  02) 


-  Sjp<W  -  “j(“p>  +  W1 


bjp  •  ("j(“p'  ■  W1  (1“p  ■  XP’ 


to  obtain  better  approximations  for  A  and  b.  for  p=l ,  ....  m  and  j=l , 

P  J  P 

m.  The  values  of  A  and  b.  so  obtained  can  then  be  used  again  in 
P  JP 


Equations  12  and  13  to  compute  still  better  values  of  A  and  b.  and  so 

P  Jr 

on  until  prescribed  tolerances  are  satisfied  (or  a  fixed  maximum  number  of 
iterations  completed). 
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Performing  the  computations  described.  Equations  10,  11,  12  and  13 
determine  vectors  b^,  k=l ,  ...,  m  where  b^  has  components  b^  for  j=l , 
...,  m.  Hence  from  Equation  2  we  have 


Vkr  =  bk  (,4) 

For  the  case  where  the  matrices  M,  C  and  K  are  symmetric,  excitation  at  a 
single  point  is  all  that  is  needed  to  determine  the  vectors  u^  and  vk; 
since  v^  is  just  some  scalar  multiple  of  uk.  Thus  we  may  take 


r  = 


(15) 


where  e^  denotes  the  coordinate  vector  with  first  component  1  and  all 
other  components  zero.  Moreover,  we  may  suppose  the  vectors  u^  normalized 
so  that  the  first  component  of  uk  has  the  value  1. 

“ik  ■ 1  061 

Then  from  Equation  14  we  have 

V  =  ?vik  '  bu 


or 


'Ik 


V 


(17) 


and  for  j  f  1 

Ujk  =  b  jk  /blk 


(18) 


and 


jk 


=  v 


lk 


jk 


(19) 


In  order  to  describe  the  process  for  the  nonsymmetric  case,  it  is 
convenient  to  modify  slightly  the  notation  used  in  the  symmetric  case. 

Let  us  denote  the  vectors  b^  by  b^O).  The  1  in  parenthesis  indicates 
that  the  vectors  b^(l)  are  determined  from  data  obtained  by  exciting 
Equation  1  at  station  1.  If  we  had  excited  Equation  1  at  station  q  with 
the  same  set  of  frequencies  as  used  at  station  1,  then  from  the  associated 
response  we  could  have  determined  vectors  b^(q)  in  exactly  the  same  way 
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(using  the  same  equations)  as  the  vectors  bk(l)  were  determined. 
However,  we  will  see  that  it  is  not  necessary  to  determine  the  vectors 
t>k(q )  for  k=l ,  m,  and  q=l ,  ....  m  completely. 


Nevertheless,  for  the  moment  suppose  that  at  the  same  time  A k  and 
b^O)  are  being  determined  iteratively  we  compute  also  b^(  2 ) ,  ...,  bk(m) 
for  k =1 ,  ....  m.  For  the  nonsymmetric  case  only  Equation  19  does  not 
hold  but  from  Equations  16,  17  and  18  we  determine  from  the  vectors  ( 1 ) 
the  vectors  uk  and  the  component  v^k  of  vk  for  k=l ,  ....  m.  Thus  we 
need  the  remaining  components  of  the  vectors  vk- 


For  the  vectors  bk(q),  for  k=l ,  ....  m  we  have  the  relation 

Vk  V  *  Vqk?  '  bk(q) 


In  this  vector  equation  there  is  only  one  unknown,  namely,  the  scalar 
v^.  From  this  it  is  clear  that  all  we  need  to  determine  the  scalar  v^k 
is  one  component  of  the  vector  b,(q)  for  k=l,  ...»  m. 


In  order  to  summarize  the  information  required  for  the  nonsymmetric 
case  let  us  rewrite  Equation  2  as 

n  Vk  eo? 

. .  /  _  \  __  r  IvfvU  /oi 


yfo.q)  =  ^  ~ ttjtt 

k=l  lu  * 


Let  Wj(w,q)  denote  the  jth  component  of  the  vector  y(u),q).  Then  in  place 
of  Equation  3  we  have 

m  _  — 

Wj  (u>,q)  =  l  (bjk(q)/(ick3-Xk)  +  bjk(q)/(iu>-Ak) )  (?2) 


For  the  nonsymmetric  case  we  need  the  values  w.(w  ,1)  and  w.(w  ,1)  for 

J  P  J  ^  P 

p=l ,  ...,  m  for  each  value  of  j  for  j=l ,  ...,  m.  We  need  also  Wj(u>p,q) 
and  Wj(u>m  +  p,q)  for  p=l ,  ...,  m  for  each  value  of  q  for  q=2,  ...,  m  but 
only  for  one  value  of  j. 
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As  in  the  symmetric  case  set 

Wjp(u),q) 

Then 


And  we  obtain  therefrom 

(iVV)Cwj(“p,q)'“jp(“p,<,)] =  (VWV1  -  "jp(Vp-q»  <25> 

The  starting  values  for  Xk  are  computed  from 

XP  ■  (CVpwj("mtp>,)  •  Vj(y')]/[wj(Vp>,!  -  ,26) 

for  any  convenient  value  of  j  for  p=l ,  ....  m.  The  starting  values  for 
b.  (1)  are  computed  from 

J  r 

bjP(1)  s  (1vV  wj(u)PJ)  (27) 

for  p=l ,  m  and  for  j=l ,  ....  m.  Starting  values  for  b.  (q)  are 

a  JP 

computed  from 

bjp(q)  =  ‘  V  (28) 

for  say  j=l  and  for  p=l ,  ....  m  and  q=Z,  m. 

Once  we  have  approximations  for  X.  and  b..(q),  the  function  w.  (w,q) 

K  Jk  JP 

can  be  evaluated  at  w  and  w  .  Improved  values  of  X.  can  then  be 

p  m  +  p  k 

computed  from  the  equation 


bjp(q)/(iu>-Xp)  +  I  (b^k(q)/(ito-Xk)  +  b  ( q )  /  ( i  uj- XR ) 
Up 


(23) 


b •  (q)/(iw  -  X  )  =  w •  (w,q)  -  w.  (w,q) 
J  P  p  J  J  P 


(24) 


X 


P 


^'m-t-p^j^m-t-p*1  \  '  ^jp^m+p’1^ 


Cwj(W,)  '  Vw'0  ' 


Wv1*  ~  Sp(V1)] 

wj(V1}  +  Vv1)] 


(29) 
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for  any  convenient  value  of  j  and  for  p=l ,  ....  m.  Improved  values  of 
b.  (1)  are  computed  from  the  equation 

J  r 

bjP(1)  =  twj(v1) '  W,,(S '  V  (30) 

for  p=l ,  m  and  j=l ,  . ..,  m.  Similarly,  improved  values  of  b.  (q) 

J  P 

are  obtained  from 

bjp(q)  =  (Wj(<yq)  -  Wjp(u)p,q))(iu)p  -  Xp)  (31  ) 

for  say  j=l ,  p=l ,  ....  m  and  q=2,  ....  m.  The  values  of  Ap  and  bjp(q) 
obtained  from  Equations  29,  30,  and  31  are  used  again  in  these  equations 
to  obtain  improved  values  and  so  on  until  tolerance  requirements  are 
satisfied  or  a  fixed  number  of  iterations  completed. 


From  the  values  of  b^(q)  obtained  from  the  above  computations  we 


now  compute 


vlk  ’  blk("/f 


for  k=l ,  . . . ,  m.  Next 


v-y^ikO) 

for  j=2,  ...,  m  and  k=l ,  ....  m.  Recall 


for  k=l ,  . . . ,  m.  Then 


(33) 

(34) 


Vqk  =  b1k(£l)/?  (35) 

for  q=2,  ...,  m  and  k=l ,  ...,  m.  It  is  clear  that  the  value  of  r  in 
Equations  32  and  35  could  be  different  for  different  values  of  the 
index  q. 
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SECTION  III 

TWO  CLOSE  CHARACTERISTIC  VALUES 

In  this  section  we  describe  a  modification  of  the  procedure  given  in 
Section  II  which  enables  us  to  handle  the  case  of  two  close  characteristic 
values.  In  order  to  describe  this  procedure  we  consider  a  complex 

valued  function  w(w)  of  the  form 

2 

w(ui)  =  E  (b./(iu>  -  A.  )  +  b./(iu  -  A.  )) 

k=l  k  (36) 


I 


The  problem  is,  as  in  Section  II,  the  determination  of  the  A^  and  b^; 
but  complicated  by  the  fact  that  and  A,,  are  close.  Hence,  for  example 
we  may  not  neglect  the  term  b2/(iw  -  A^)  when  w  is  close  to  Im[A^]. 

First,  we  want  to  make  clear  the  relation  of  the  function  w(m)  of 

Equation  36  to  the  frequency  response  function  and  the  iteration  process 

of  Section  II.  Thus  w(ui)  denotes  the  value  of  the  jth  component  of  the 

frequency  response  function,  that  is,  the  function  w.(cj)  of  Equation  3 

J 

with  all  the  terms  for  k  >  2  either  ignored  as  a  first  approximation  in 
the  iteration  process  or  subtracted  off  for  the  succeeding  steps  of  the 
iteration  process.  The  subscript  j  has  been  left  off  for  notational 
simp! icity. 


Set 


w(u)  =  w(u>)  -  b^/(iu>  -  A-j )  -  b^/ ( ico  -  A^) 


(37) 


that  is 


w ( ui )  =  b1/(iu>  -  A-| )  +  b2/(iw  -  A^) 


(38) 


We  suppose  w^.  =  w(uk)  is  known  for  four  values  of  u>.  close  to  Im[A^]  and 
to  ImU^].  Using  this  data,  Equation  38  becomes  a  system  of  four 
equations  for  the  four  unknowns  b-|  ,  b^,  A^  and  A2>  The  unknowns  b^  and 
b2  are  easily  eliminated  from  this  system.  In  doing  so  one  obtains 
two  equations  for  A-j  and  A2>  We  will  now  give  some  of  the  details. 


V 
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Solve  Equation  38  for  b-j 

b,  =  w.(iu).  -  X.)  -  b9(iw.  -  X,  )/(ico .  -  X,)  (39) 

1  J  J  I  ^  j  1  J  4 

then  subtract  Equation  39  for  j  t  4  from  Equation  39  with  j  =  4  and 
after  some  rearranging  and  simplifying  obtain 

ib2(X]  -  X2)  =  [i(wjwi  -  u4w4)  -  X-j (Wj  -  w4)]P(2,  j,  4 )/{u.  -  w4)  {40) 
where 

P(k,  j.  O  8  ^  *  i(“j  +  ‘  (41) 

for  l  =  4  and  j  =  1,  2,  3;  that  is,  a  system  of  three  equations  for  the 
unknowns  b^,  X^  and  X2<  Clearly,  b^  is  readily  eliminated  from 
Equation  40  leaving  us  with  a  system  of  two  equations  which  we  can  write 
in  the  form 


f(w3  -  w4)P(2,  3,  4) 


(w.  -  w4)P(2,  j,  4)1 


(42) 


(W3w3  '  W4W4^P^2,  3>  4)  (^jwj  ~  tjLi4w4 ) p  ( 2 »  J>  4) 


w4  '  w3 


03.  -  03  . 

4  j 


for  j  =  1 ,  2 . 

Observe  that  X^  occurs  linearly  in  this  system  of  equations.  Hence 

Xj  could  be  eliminated  leaving  a  single  equation  of  the  fourth  degree  in 

X2 .  This  resulting  polynomial  equation  could  be  solved  for  X2  and,  of 

course,  we  would  have  to  choose  the  correct  root.  Thus  taking  w.  =  w(oj.) 

J  J 

as  a  first  approximation  we  could  determine  X2  as  described,  compute  X^ 

.  from  Equation  42,  b2  from  Equation  40  and  b^  from  Equation  39.  In  this 

way  one  obtains  starting  values  for  b^ ,  b2>  X ^  and  A,,.  These  initial 


V 
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values  for  b^  ,  b^,  and  A^  can  be  used  in  Equation  37  to  compute  better 
values  of  w.  for  j  =  1,  2,  3,  4  and  the  procedure  just  described  for 

J 

calculating  A£»  A-j ,  b^  and  b^  repeated. 

Since  solving  a  fourth  degree  polynomial  for  A^  and  selecting  the 
appropriate  root  is  not  an  attractive  procedure,  we  sought  an  alternative. 
Let  us  make  two  observations.  First,  if  A£  is  known  then  Equation  42 
uniquely  determines  A-j .  Secondly,  if  we  had  first  eliminated  b£  and  then 
b.|  we  would  have  obtained  Equation  41  with  the  roles  of  A^  and  A^ 
interchanged;  so  if  we  know  A^  then  A^  is  uniquely  determined.  All  the 
equations  needed  for  treating  the  case  of  two  close  characteristic  values 
are  now  available. 


We  do  not  feel  that  we  have  sufficient  numerical  experience  to  give 
a  rigidly  fixed  recipe.  However,  of  the  various  alternatives  which  we 
tried  the  procedure  which  we  describe  next  seemed  best.  We  assumed  an 
initial  approximation  to  Ap.  From  Equation  41  we  compute 

2 

P  ( 2 ,  j,  4)  =  Ap  ~  i  ( ui j  +  u>4 ) A 2  ~  w j w4  ,  4 

for  j  =  1  and  3.  Thc-n  from  Equation  42  we  have 


/V 


(w3  -  w4)P(2,  3,  4)  (W]  -  w4)P(2,  1,  4) 


J4  '  u'3 


w4  - 


(44) 


-  4w4)r(2,  3,  4)  (.^  -  ,'4w4)P(2,  1,  4) 


''4  "  ul3 


w4  - 


which  is  easily  solved  for  A^ .  Now  that  A^  and  Ap  are  known  we  use 
Equation  40  to  compute  bp,  Equation  39  to  compute  b-j  and  Equation  37  to 
obtain  improved  values  of  wn  for  j  =  1,  2,  3,  and  4. 
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Next  we  compute 

P(l,  j.  1)  = 

for  j  =  2,  4  and  then  solve 


1  ((Jjj  +  U)-|  )  A-|  - 


(45) 


(w2  -  w^Pd,  2,  1)  (w4  -  Wl)P(l,  4,  iy 


~ '  CO  2 


w-|  -  w4 


(46) 


(^2^2  ~  w-|wi)p(l»  2»  1)  (<^4W4  -  w.|Wj)P(l,  4,  1) 


-  tog 


“1  ‘  “4 


for  A2>  Now  th"*t  a  new  value  of  A2  is  available  we  use 

A  A  A  A 

^  i  b  1  ( X  2  -  ^  )  =  [i(w.ujj  -  w-|tOi )  -  A2(Wj  -  w] )]  P(l,  j,  l)/(wj  -  u)-, )  (47) 

with  either  j  =  2  or  j  =  4  to  compute  b-|  and 

•  * 

j  b?  =  Wj ( i uk  -  \2)  -  b^iu-^  -  A2)/(i^  -  A-j )  (48) 

again  with  either  j  =  2  or  j  =  4. 

In  the  procedure  described  above  we  have  assumed 

0  <  Im[A-|]  <  Im[A2] 
and 


cj-|  <  oj^  <  0)^  <  ^4  . 

It  also  seemed  better  to  perform  two  successive  updates  of  the  close 
characteristic  values  and  their  associated  coefficients  for  every  update 
of  the  well  separated  characteristic  values  and  their  associated 
coefficients,  at  least  in  the  later  stages  of  the  iteration  process. 
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SECTION  IV 

THE  FREQUENCY  RESPONSE  FUNCTION  BY  NEWTON'S  METHOD 

In  Section  II  we  described  the  determination  of  the  frequency 
response  function  from  experimental  data  by  an  iteration  process.  In 
this  section  we  describe  the  determination  of  the  frequency  response 
function  by  Newton's  (iteration)  method.  The  method  of  Section  II  and 
the  Newton  method  are  not  too  different  from  one  another. 

For  simplicity  we  leave  off  the  subscript  denoting  the  component  of 
the  frequency  response  function  (Equation  3).  The  problem  is,  as  in 
Section  II,  the  determination  of  the  unknowns  b^  and  in  the  function. 

m  . 

w(u)  =  Z  [b^/  ( ico-Xk )  +  bk/(ito-Ak)]  (49) 

from  values  w.  =  w(w.)  for  j=l ,  ...,  n=2m. 

J  3 

Set 

m 

W(ui,  b-j ,  A-| ,  ...,  b^,  >>m)  =  ^  bk/(iw-Ak)  (50) 


W(w,  b.|  ,  ,  ...,  b^,  A  )  is  a  function  of  the  variables  to  ,  b^  ,  A^ . 

b  ,  A  ,  Suppose  we  know 
m  m 

Wj  =  W(w.,  br  Ar  ...,  bm,  Am)  (51) 


for  n  values  of  j  and  we  have  some  initial  "guess"  b^°,  A^°,  ...,  bm°, 

A  0  for  the  values  bn ,  A,,  ...,  b  ,  A  .  Then  from  Equation  50  we  can 
m  11mm  M 

compute 


“d  ■ 


b;. 


1° 
A1  ’ 


Am> 


(52) 


for  J=1 ,  ...,  2n  and  from  the  system  of  equations 


3W  ..  ,  9W_  hi  + 

3b1  dbl  9A]  dAl 


+  fjj-  db  +  It-  dA  =  W.-WY 
9bm  m  m  ■  J  J 


(53) 
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for  j=l . n  we  can  determine  the  increments  db, ,  dX.  ,  ....  db  ,  dA 

1  1  m  n 

for  improving  the  approximations  b^°,  X^,  ....  b°,  X®  .  Thus 

b°  (new)  =  b°  +  dbfe 


X®  (new)  =  X°  +  dXk 

In  the  system  of  equations  defined  by  Equation  53 

3W/3bk  =  l/(ia>j  -  A°) 
and 

3W/3Xk  =  b®/(iu)j  -  X®)2 

In  the  practical  situation  we  do  not  have  the  function  values  W. 

J 

but  have  instead  the  function  values  w..  Hence  we  replace  W.  in  the 

J  .1 

system  Equation  53  by 

m  _ 

W,  =  w,  -  E  b P / ( i to  .  -  X?) 


J  J 


k  j  Aky 


The  initial  guess  for  the  bk  and  Xk  is  obtained  in  the  same  way  as 
in  Section  II.  The  main  difficulty  with  this  process  is  solving  the 
n  by  n  system  of  complex  equations.  Equation  53,  for  each  iteration  step. 
At  this  time  we  believe  the  method  of  Section  II  will  involve  fewer 
computations  even  though  it  may  require  more  iterations.  The  Newton 
process  given  here  is  for  the  case  of  well  separated  characteristic 
values  only. 

Here  now  we  re-examine  briefly  the  procedure  presented  in  Section  II. 
We  can  simply  the  notation  and  still  illustrate  the  main  features  if  we 
consider,  as  above  Equation  49,  a  single  complex  valued  function.  We 
suppose  the  value  w(to)  given  for  2m  values  of  to.  Thus  we  take 


w(to.)  =  w. 

J  J 

W  ( (0  .  .  )  =  w  .  . 
m+j '  m+j 
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for  j  =  1 ,  ....  m.  We  suppose 

<  Im  [Xj3  <  Vi  (59) 

and  that  the  difference  gj  ..  -  g>.  is  small.  The  problem  is  the  same  as 

m+j  j 

above,  namely,  solving  the  resulting  system  of  equations  for  the  m 
coefficients  b^  and  the  m  characteristic  values  A^. 

For  the  2m  values  of  ui  the  associated  system  of  equations.  Equation  49, 
is  nonlinear  in  the  A^.  The  first  step  in  the  solution  process  is  to 
rewrite  this  system  of  equations.  Set 


W.  =  bj/(ioK-Xj)  +  £  [b^/OaK-A^)  +  b^Ou^-A^)]  (60) 

k/j 


Wm+j  is  defined  similarly. 


The  system  of  Equation  49  now  appears  as 
w.  =  b./(iu»rA.)  ♦  Mj 

wm+j  =  bj^1wm+j'Xj^  +  Wm+j 


(61) 


for  j=l ,  ...»  m.  For  a  fixed  value  of  j  the  two  equations  in  the  system. 

Equation  61,  can  be  solved  for  A.  and  b..  We  have 

J  0 

Xj  =  'tVjVT  "j"j  •  <VlVi  -  “jV]/tVi  '  "j  •  ('Vj  ■  “j’3 

(62) 

bj  *  ("j  '  Wj)(1“j  '  ‘j1 


for  j  =  1,  ...,  m.  Now  then  let  us  regard  the  A.  and  b.  as  variables. 

J  J 

If  we  substitute  in  the  right  hand  side  of  Equation  62  values  A.  ,  b. 
for  j  =  1 ,  ....  m  we  will  obtain  say  values  A ^ ,  b^ ,  for  j  =  1,  ...,  m. 
Thus  the  system  of  Equation  62  transforms  a  set  of  2m  complex  numbers 
into  a  set  of  2m  complex  numbers.  The  desired  solution  is  a  set  of  2m 
complex  numbers  which  transform  into  the  same  set  of  complex  numbers 
respectively. 
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In  Section  II  we  take  b,  =  . . .  =  b  =0  initially.  For  this 

i  m 

particular  assignment  of  the  coefficients  bk  it  is  not  necessary  to 
specify  the  Afc.  Then  for  j  =1  we  use  Equation  62  to  compute  and  b1 . 
We  could  have  computed  A.  and  b.  for  j  =  2,  ....  m  also  for  b,  =  ...  = 

J  J  * 

b  =0.  Instead,  we  use  A,,  b,  and  b0  =  ...  =  b  =  0  to  compute  A,  and 
m  l  l  i  m  Z 

b^  and  so  on. 

We  would  like  to  be  able  to  show  that  the  transformation  defined 
by  Equation  62  is  a  contraction  mapping.  Then  it  would  follow  that  the 
iteration  process  of  Section  II  converges  to  a  unique  fixed  point. 

As  of  this  time,  however,  we  have  not  been  able  to  determine  a  region 
of  convergence. 
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SECTION  V 

RESULTS  AND  CONCLUSIONS 

In  this  section  we  describe  and  discuss  briefly  the  numerical 
experiments  performed  to  test  the  procedures  presented  in  Sections  II  and 
III.  From  the  discussion  there  and  Appendix  A  the  problem  of  interest 
reduces  to  the  determination  of  parameters  b^  and  A^  which  define  a 
complex  valued  function  w(ta)  of  the  form 

m  _ 

w(oj)  =  2^  [bk/(ito  -  A^)  +  b^/dw  -  Ak)] 

Our  first  objective  is  to  exhibit  the  influence  of  a  coefficient  b^, 
in  our  case  b^,  on  the  function  w(ui).  The  characteristic  values  Afc  and 
coefficients  b^  used  for  defining  functions  w(u>)  are  listed  in  Tables 
1A-4A.  The  values  of  A^  and  b^  are  the  same  in  all  these  tables  except 
for  b^.  Since  Ak  and  b^  are  complex  quantities  it  takes  a  pair  of  real 
numbers  to  specify  them.  The  functions  w(a>)  resulting  from  simple 
modifications  of  the  coefficient  t>4  are  shown  graphically  in  Figures  7a-10a. 
Figures  7b-10b  are  the  graphs  of  |w(u>)|  for  the  same  modifications  of 
the  coefficient  b^.  In  order  to  facilitate  comparison  these  graphs  are 
shown  collectively  in  Figures  6a  and  6b  respectively. 

One  observes,  Figure  6a,  as  the  coefficient  b^  changes  from  quadrant 

to  quadrant  there  is  a  corresponding  change  in  one  loop  of  the  associated 

graph.  The  remaining  portion  of  the  graph  is  not  visibly  changed.  Thus 

it  is  clear  that  the  term  b^/d'a-A^)  is  the  principal  contributor  to  the 

value  of  the  function  w(w)  when  is  in  the  vicinity  of  I  [A.].  This 

J  m  k 

feature  is  illustrated  further  in  Figures  7a-10a  where  the  x's  are  the 
values  of  the  term  b^/d’w-A^)  alone  as  w  varies  from  4.3  to  5.7  at  steps 
of  0.05. 

In  Table  1A,  for  example,  twelve  frequencies  (in  two  columns)  are 
listed.  The  corresponding  function  values  of  w(w)  are  listed  immediately 
thereafter.  These  function  values  are  indicated  on  the  corresponding 
graphs.  Figures  7a  and  7b,  by  the  "boxes".  From  these  twelve  function 
values,  using  the  method  of  Section  II,  the  characteristic  values  A^  and 
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coefficients  bk  are  determined.  The  characteristic  values  and  coefficient 
computed  from  the  given  function  values  are  listed  in  Table  IB.  The 
"radius  of  differences"  is  the  square  root  of  the  sum  of  squares  of  the 
real  and  imaginary  parts  of  the  differences  of  the  last  two  successive 
calculations  of  the  characteristic  values. 

Table  1A  also  gives  "truncated"  values  of  the  function  w(u>)  at  the 
listed  frequencies.  Table  1C  gives  the  characteristic  values  and 
coefficients  obtained  from  computations  using  these  truncated  values. 

Thus  Table  1C  shows  how  much  the  lack  of  precision  in  the  experimental 
data  affects  the  computations.  Tables  2A-4C  show  the  results  as  the 
coefficient  b^  is  changed.  As  one  would  expect  the  method  of  Section  II 
is  not  sensitive  to  changes  in  the  coefficient  b^. 

Our  second  objective  is  to  test  the  procedures  described  in  Section  III 
for  two  close  characteristic  values.  The  data  and  results  for  this  set 
of  problems  is  given  in  the  Tables  5A-9D.  The  coefficient  b^  and 
characteristic  values  A^,  except  for  Ag,  are  the  same  throughout  the  set 
of  problems.  The  characteristic  value  Ag  gets  close  to  Ag  by  letting 
Im[Ag]  take  on  the  values  4.5,  4.3,  4.2,  4.1  and  4.05  successively. 

In  order  to  see  the  effect  of  Ag  tending  to  Ag  the  graphs.  Figures  12a- 
15a  and  12b-15b  are  collected  in  Figures  11a  and  lib  respectively. 

One  observes,  in  Figure  11a,  that  as  Ag  tends  to  Ag  the  "circle" 
corresponding  to  Ag  deteriorates,  diminishes  and  finally  disappears.  At 
the  same  time  the  "circle"  corresponding  to  Ag  distorts  and  expands  and 
eventually  becomes  nearly  circular  again  when  Ag  is  practically  the 
same  as  Ag.  Similarly,  the  valley  between  the  peaks  associated  with  Ag 
and  Ag  (Figure  lib)  disappears  as  Ag  tends  to  Ag  and  the  two  peaks  become 
a  single  high  peak  when  Ag  is  practically  the  same  as  Ag. 

It  is  clear  that  the  behavior  just  described  depends  upon  the 
coefficients  bg  and  bg.  To  illustrate  the  role  of  bg  and  bg  envision 
what  takes  place  as  Ag  tends  to  Ag  for  the  case  bg  =  -bg. 
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In  Tables  5B-9B  we  give  the  results  using  the  procedures  described 
in  Sections  II  and  III  and  the  exact  function  values.  In  Tables  5C-9C 
we  give  the  results  when  the  truncated  function  values  were  used.  We 
observe  that  with  accurate  data  all  coefficients  b^  are  computed  well. 

With  inaccurate  data,  the  difference  between  the  actual  and  computed 
coefficients  gets  large  particularly  for  the  differences  corresponding 
to  A,.  and  Ag  respectively  as  Ag  tends  to  Ag. 

Using  the  computed  characteristic  values  and  coefficients  from 
Table  9C  we  computed  the  function  values  at  the  given  set  of  frequencies. 
These  computed  function  values  are  listed  in  Table  9D.  Observe  that  the 
computed  function  values  in  Table  9D  and  the  truncated  function  values 
in  Table  9C  agree  very  well.  Observe  also  that  graphically,  in  Figure  17, 
the  function  values  based  on  the  computed  coefficients  from  Table  9C  are 
indistinguishable  from  function  values  based  on  the  exact  characteristic 
values  and  coefficient. 

We  believe  the  numerical  experiments  performed  in  Reference  1  and 
this  report  show  that  the  procedure,  described  in  Section  II  of 
Reference  1  and  Sections  II  and  III  of  this  report,  is  a  highly  accurate 
method  for  determining  the  complex  characteristic  values  and  associated 
complex  characteristic  vectors  of  a  light,  viscously  damped  linear  system 
from  the  responses  to  sinusoidal  excitations.  The  method  as  presented 
in  Reference  1  was  capable  of  handling  well  separated  characteristic 
values  and  characteristic  values  of  multiplicity  great  than  1. 

In  Reference  1  we  tentatively  assumed  that  the  characteristic  values 
could  be  regarded  either  as  well  separated  or  identical.  However,  we 
were  concerned  that  this  assumption  might  not  be  realistic.  Accordingly 
we  extended  the  method.  Section  III,  so  that  it  could  handle  the  case  of 
two  characteristic  values  which  were  too  close  for  the  well  separated 
procedure  yet  not  close  enough  to  be  regarded  as  identical. 


V 
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Further  numerical  experiments  should  be  performed  to  determine  more 
precisely  when  two  characteristic  values  should  no  longer  be  regarded 
as  well  separated  but  rather  as  close  and  when  no  longer  as  close  but 
identical.  We  need  also  to  perform  numerical  experiments  to  determine 
how  many  characteristic  values  the  method  can  handle.  It  is  also  of 
interest  to  determine  the  degree  of  damping  at  which  the  method  is  no 
longer  reliable. 


28 


AFWAL - TR-80-31 36 


APPENDIX  A 

LINEAR  SYSTEMS  OF  ORDINARY  DIFFERENTIAL  EQUATIONS 

Methods  for  determining  vibration  parameters  are  based  on  facts  and 
relations  from  the  theory  of  linear  systems  of  ordinary  differential 
equations  with  constant  coefficients.  These  facts  and  relations  are 
well  known,  at  least  to  the  expert.  In  this  appendix,  it  is  our  objective 
to  develop  in  a  brief,  yet  complete  fashion,  these  relations  so  that 
they  are  readily  accessible  and  acceptable. 

Consider  a  second  order  system  of  differential  equations 

Mx  +  Cx  +  Kx  =  f  (£1 ) 

Here  M,  C  and  K  denote  real  matrices  of  order  m  and  x  =  x ( t )  and  f  =  f(t) 
are  m-dimen$ional  vector  functions  of  t.  The  vector  functions  x(t)  and 
f(t)  may  have  complex  components. 

The  associated  homogeneous  equation 

Mx  +  Cx  +  Kx  =  0  (A2) 

has  nontrivial  solutions  of  the  form  u  exp(At)  if  and  only  if  A  satisfied 
the  characteristic  equation 

det[MA2  +  CX  +  K]  =  0  (A3) 

and  the  vector  u  satisfies  the  condition 

[MX2  +  CX  +  K]  u  =  0  (A4) 

Similarly,  if  a  vector  v  satisfies  the  condition 

[mV  +  CTX  +  KT]  V  =  0 

then  v  exp(Xt)  satisfies  the  transposed  homogeneous  equation 

MTx  +  CTx  +  KTx  =  0 


(A5) 


(A6) 
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Throughout  this  appendix  the  symbol  v  is  reserved  for  quantities 
associated  with  the  transposed  equation.  We  will  refer  to  values  of  A 
which  satisfy  the  characteristic  equation.  Equation  A3,  as  characteristic 
values  and  the  vectors  u  and  v  as  corresponding  (or  associated)  charac¬ 
teristic  vectors. 

The  characteristic  equation  is  a  polynomial  equation  in  A.  Since 

the  coefficients  are  real,  the  complex  conjugate  of  any  complex  root  is 

also  a  root.  We  suppose  Equation  A3  is  of  degree  n  =  2m  and  that  all  the 

roots  are  complex  with  negative  real  part.  We  suppose  also  that  if  A  is 

2 

a  root  of  multiplicity  p  then  the  matrix  [MX  +  CA  +  K]  is  of  rank  m-p. 

For  convenience  we  index  the  roots  of  Equation  A3  so  that  0  <  Im[A^]  < 
IM[\2]  -  ...  _  Im[\  ]  and  "a^  =  Afc ,  for  k  =  1  ,  . . . ,  m.  Let  uk  and  vfc 
denote  characteristic  vectors  corresponding  to  A^  for  any  value  of  k  <_  m. 
That  is,  u^  and  vk  satisfy  the  conditions  expressed  by  Equations  A4  and 
A5  respectively.  Then  u^+k  =  uk  and  vm+k  =  vk  are  characteristic 
vectors  corresponding  to  Am+k  =  Ak. 

Orthogonality  conditions  are  well  known  and  convenient  normalizing 
relations  are  readily  determined  for  first  order  systems  of  differential 
equations.  For  these  reasons  we  consider  the  first  order  system 

x  i  ”  x  ^  ~  0 

(A7) 

Mx9  +  Kx,  +  Cx0  =  f 

C  I  c 

which  i>.  a  system  equivalent  to  Equation  A1 .  This  system  may  be  written 
in  block  matrix  form  as 


I  0 


i_0  KJ  [  x 
and  symbol ically  as 


0  -I 


Ay  +  By  =  g 
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One  can  readily  verify  u^expU^t) ,  where 


Uk  =  X 


(AlO) 


satisfies  the  homogeneous  equation  Ay  +  By  =  0.  Similarly,  vk  exp(A^t) 

W'W\]  Ml 


(All) 


satisfies  the  transposed  homogeneous  equation  ATy  +  BTy  =  0.  The  vectors 
uk  and  vk  are  characteristic  vectors  satisfying  Equations  A4  and  A5 
respectively. 

From  the  equations  X^A  uk  +  B  uk  =  0  and  \  .vTa  +  vTb  =  0  it  follows 

^  **  J  J  J 

that 

V/A  \  +  B  Qk  =  0 

X  vTA  u,  +  vT  B  u.  *  0 

j  j  K  J  K 


On  subtracting  the  second  of  these  two  equations  from  the  first  one 

finds  for  j  ^  k,  if  X.  f  X.  that  the  characteristic  vectors  u.  and  v. 

J  K  *  J 

satisfy  the  orthogonality  conditions 


vj  A  u^  =  0  and  vT  B  uk  =  0 


(A12) 


For  a  characteristic  value  Xk  of  multiplicity  p  >  1,  the  characteristic 

vectors  can  be  determined  in  pairs  uk,  vk,  uk+1 ,  vk+] . uk+p,  vk+p 

so  that  the  orthogonality  conditions  expressed  by  Equations  A12  hold 
whenever  j  f  k. 


For  reasons  which  will  be  apparent  shortly,  we  suppose  the  vectors 
uk  and  vk  normalized  so  that 


*k  A  Dk  =  1 


(A13) 


_  iii. 
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Equation  All  to  replace  the  row  vectors  vj^.  Lastly,  use  the  right  hand 
side  of  Equation  A8  to  replace  g.  Thus  the  integrand  of  Equation  A18 
becomes 


expXiT  0 

( - 1 / X i ) vT  ^  K  v{ 

0 

0  expX  t 

_  n  - 

f(x) 

m 

Performing  the  indicated  matrix- vector  multiplications  this  integrand  is 
rewritten  as  the  vector 

(exp>.-|2)  v ^ Tf (t ) 

(expXn2)  vjf(2) 


Multiplying  this  vector  by  the  matrix  UAe(t)  Equation  A18  can  be 
rewritten  as  the  two  vector  equations 

n  f1"  T 

x-j  ( t )  =  I  uk  exp(Xkt)  I  exp(-XkT)vk  f(x)dT  (A19) 

k=1  JQ 


n 

x-(t)  =  Z 
1  k=l 


Xkuk  exp(Xkt)  I  exp(-XkT)vk  f(x)dT 


(A20) 


It  follows  that  Xl(t)  and  x2(t)  must  satisfy  Equation  A7.  From 
Equation  A7  x^(t)  =  XjU).  If  we  compute  Xj(t)  from  Equation  A19,  we 
find  that 


x. (t)  =  x-(t)  +  Z  u.vT  f(t) 
1  £  k=l  K  K 
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and  it  follows,  since  f(t)  is  arbitrary,  that 

I  u.  vj  =  0  (A21 ) 

k=l  K  k 

In  a  similar  way,  we  find,  since  x^t)  must  satisfy  the  condition 
Mx  +  Cx  +  Kx  =  f,  that 

M  •  l  AkukvJ  =  I  (A22) 

k-1 


Using  Equation  All  we  can  write  Equation  A16  as 


-kt  0 

Vl/Xl  •••  vn/  n 

0  I_ 

1 

< 

. — • 

< 

1 - 

( A2  3 ) 


then 


VTA  U 


Tm 

Vl/Xl 


v]/X 
n  n 


-K  0 


0  M 


u,  . ..  u. 


Lx,u 


r  i 


••  Xu 

n  n 


(A24) 


In  the  same  way  we  have 


vtb  u  = 


r 


v^/X 

n  n 


0  K 


U  C. 


u,  . . .  u 


Lx,u, 


X  u„  . 
n  n  J 


=  -  A 


(A25) 


Equations  A22,  A24  and  A25  show  that  the  characteristic  values  X1 ,  ..., 
Xm  and  appropriately  normalized  characteristic  vectors  u^ ,  . . . ,  u^  and 
V| ,  . . . ,  vm  determine  the  matrices  M,  C  and  K.  Actually,  Equations  A22 
and  A25  determine  M,  C  and  K.  The  Equation  A25  is  a  matrix  equation  of 
order  n  =  2m.  Other  relations  are  available  for  determining  the 
matrices  M,  C  and  K. 
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It  follows  from  VTAU  *  I  that  AUVT  <*  I.  Using  this  last  equality 
one  can  obtain  the  results  expressed  by  Equations  A21  and  A22.  With 
Equation  A21  and  Equation  A22  established  in  this  second  Independent 
manner,  it  is  a  simple  task  to  verify  that  the  time  derivative  of 
Equation  A19  is  Equation  A20.  It  is  also  readily  verified  that  ( t )  as 
given  by  Equation  A19  satisfies  Equation  A1 . 


From  Laplace  transform  theory  we  have  for  zero  initial  conditions 

L{x(t) }  =  [Ms2  +  Cs  +  K]"1  L{f (t)}  (A26) 

We  also  have 

ft 

L{J0  Fi(t-T)F2(T)dT}  -  f,(s)  •  f2(s)  (A27) 

where  f^s)  =  L{F-j  ( t ) >  and  f2(s)  =  L{F2(t)}.  Applying  the  formula 
expressed  by  Equation  A27  to  Equation  A19  we  obtain 


L{xl(t)}  »  E  ^r-Uf(t)} 
k=l  k 

From  Equations  A26  and  A28  we  see  that  the  transfer  function 

9  .1  ** 

[Ms£  +  Cs  +  K]  '  «  E  -K~  «  H ( s ) 
k*l  s"Ak 


(A28) 


(A29) 


say. 


Next  we  want  to  determine  some  particular  solutions  of  Equation  A1 
for  some  simple  vector  functions  f(t).  Let  ej  denote  the  vector  whose 
jt.h  component  Is  1  and  all  other  components  zero.  Also  let  h(t)  denote 
the  scalar  function  satisfying  the  conditions. 


h(t)  =0  for  t  <  0 
=  1  for  t  >  1 
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Set  f(t)  =  e j h ( t )  and  denote  the  response  to  ejh(t)  by  H ^ ( t ) .  Replacing 
f(t)  in  Equation  A1  9  by  e.h(r)  one  immediately  obtains  for  the  step 

J 

response  H.(t) 

J 


"  Jkvkej 


!  Vk*j 


Hj ( t)  =  I  -V-J-  exp(X.t)  -  E  ~ 
k=l  Ak  K  k=l  A 


(A30) 


The  impulse  response  I.(t)  is  the  derivative  of  the  step  response,  that 

J 


Ij<‘>  -df  ■  j,  Vh  <V>  (A31 ) 


Let  fj(T)  denote  the  jth  component  of  the  vector  forcing  function 
f(t).  Then  f(t)  can  be  written  as 

f(t)  -  e,f,(t)  *  *  emfm(t). 

Replacing  f(x)  in  Equation  A19  with  this  expression,  we  obtain 

m  t  n 

x,(t)  =  I  f  E  expx.  (t-r)u.  v,  e.f.(x)d 
1  j=l  k=l  k  k  k  J  J 


From  Equation  A31  this  formula  for  the  response  can  be  written  as 


m  t 


<l(t)  -  I  f  Ij(t-T)  fj(T] 
J-1  •'n 


(A32) 


Next  we  take  f(t)  =  r  exp(iust)  where  r  denotes  a  real  constant  vector 
and  obtain  from  Equation  AT 9 

n  uifVmr  n  UlkVtr 

x,(t)  ‘  k=i  e>"’(1“t>  ‘  k=i  ^ exp(>kt>  (A33) 


For  complex  Xk  with  negative  real  part,  the  second  summation  in  Eq.  A33 
goes  to  zero  as  t  becomes  large.  Hence  ^ 


"  Vkr 


y  exp(iu)t)  =  E  ■  — t  -  exp(iwt) 
k=l  k 


(A34) 


is  the  steady  state  response  to  the  harmonic  excitation  r  exp(i<ut). 
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The  steady  state  response  to  r  exp(iu)t)  is  readily  obtained  from 
the  steady  state  response  to  the  sinusoidal  excitation  r  sin  ut.  One 
knows,  or  can  readily  show,  that 

[MD2  +  CO  +  K]  {In[y  exp(iwt)]}  =  Im  [r  exploit)]  =  r  sin  u>t  (A35) 
Now 

Im[y  exp(ia)t)]  =  Re  [y]  sin  tot  +  Im[y]  cos  u>t  (A36) 

and  if  x(t)  denotes  the  steady  state  response  to  the  r  sin  cot  the  system 
of  equations 

Re[y]  sin  tot^  +  Im[y]  cos  u>tj  =  x(t-|) 

\  A3/ ) 

Re[y]  sin  wt2  +  Im[y]  cos  wt2  *  x(t2) 

is  readily  solved  for  Re[y]  and  Im[y]  from  the  recorded  values  of  x(t) 
at  t1  and  t2,  provided  sin  tot1  cos  tot2  -  cos  wt^  sin  ut2  *  0. 


From  Equation  A34 


n 

y  *  l 

k=l 


i  a.'-  X 


n 

[  z 

k=l 


■]r 


(A38) 


The  matrix  on  the  right  hand  side  of  this  equation  is  called  the 
frequency  response  function.  Observe  that  it  differs  from  the  transfer 
function,  Equation  A29,  only  in  the  term  i to  in  the  denominator. 


From  Equation  A38  we  have 

n  u.vJr  m  u.vJr 

y  =  l  +  Z  -- 

k=l  -iw-X  k=l  - i oj- X 

and  from  this  equation  it  is  clear  that 

n  u.vJr 
y  ■  z  -f 


k=l  'iw'Xk 


(A39) 
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Consider  next  the  undamped  homogeneous  system  of  equations 

Hx  +  Kx  =  0  (A4°) 

associated  with  Equation  A1 .  We  now  suppose  that  M  and  K  are  symmetric 

and  that  the  positive  definiteness  condition  satisfied.  Then  there  are 

values  0  <  <  . . .  £  and  corresponding  real  rectors  u^ , 

u  satisfying  the  condition 
m  J 

1-4  M  +  K]  uk  =  0  (A41 ) 

It  is  readily  shown  that  the  vectors  uk  satisfy  the  orthogonality 
conditions 

"'T  ~  ~T 

UjMuk  =  0  and  u.  K  =  0  (A42) 

We  will  not  impose  any  particular  normalization  on  the  vector  uk  at 
this  time. 


Set 

A  rt  A 

-irk  exp ( i ' >k t )  =  [MD^  +  CD  +  K]  (uk  exp(icokt)} 
and 

A  a  A 

irk  exp  (-icokt)  =  [MD^  +  CD  +  K]  (uk  exp(-iu>kt)} 
On  adding  these  two  equations  one  obtains 

9  ~  ~ 

(MD  +  CD  +  K]  ( u k  cosu>kt}  =  rk  sin^t 
This  equation  simplifies  to 

'wk  C  uk  =  rk 

for  k  =  1,  ....  m.  This  system  in  matrix  form  is 


(A43) 


C[U,  ...  um] 


■U>1  0 


0  -co 


h  •••  --J 


(A44) 


V 
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Equation  A44  can  be  solved  for  C  and  we  have 


C  “  Crl  rm] 


-1/c^  0 


■1/^ 


C«1 


(A45) 


Set 


U  =  [Ul  ...  um] 

Because  of  the  orthogonality  conditions 


(A46) 


AT  A 

U  MU  = 


al 
0  a, 


m 


(A47) 


A  diagonal  matrix  with  diagonal  entries  a. . a_  which  must  be 

I  Ml 

determined.  Similarly 

2,a,  0  1 


aT  a 

u'ku  = 


u 


ri 


0  u,2a 
m  m 


Next  for  some  value  of  w  t  we  have 

[MD2  +  CD  +  K]  y  exp(iwt)  =  r  exp  (iut) 

From  this  equation  one  obtains,  since  we  may  suppose  C  Is  known 


[-<o  M  +  K]y  =  r  -  iuCy 


Then 


Set 


UT[-0J2  M  +  K]UU"'y  =  U'(r  -  iwCy)  =  b  say 


U'V  - 


lzn 


(A48) 


(A49) 
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If  we  examine  Equation  A52  and  A54  together  we  are  able  to  observe 
an  important  relation  between  the  transfer  function  (or  better  the 
frequency  response  function)  of  Equation  1  and  the  natural  modes  and 
natural  frequencies  of  the  associated  undamped  system.  Suppose  w  =  u>k  is 
a  characteristic  value  of  Re[H(iu>)]  and  r  is  a  corresponding  characteristic 
vector,  that  is 

det[Re(Hitt>k)]  =  0 
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and 

Re[y]  =  Re[H(iwk)]r  =  0 

Now  the  vector  y  determined  from  Equation  A54  by  r  and  wk  must  satisfy 
Equation  A52  and,  in  particular,  since  y  =  i  Im[y]  it  follows  that 

[ M  +  K]  y/i  =  0 

Hence  y/i  and  y  are  characteristic  vectors  associated  with  the  character¬ 
istic  value 

Conversely,  if  u  =  a>k  and  a  real  vector  which  we  denote  by  Im[y] 
satisfy  the  condition 

[-wk  M  +  K]  Im[y]  =  0, 

then  from  Equation  A52  we  infer  Re[y]  =  0  and  r  =  Im[y].  Since  these 
quantities  must  satisfy  Equation  A54,  it  follows  in  particular  that  r 
satisfies  the  homogeneous  equation  Re[H(iu>k)]  r  =  0.  Hence  satisfies 
det[Re[H(i<jok)]  =  0. 
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APPENDIX  B 

INTERPOLATION  WITH  EXPONENTIAL  FUNCTIONS 


In  this  section  we  give  the  mathematical  background 
some  procedures  for  determining  the  "amplitudes"  ak  and 
frequencies"  Xk  of  a  function  x(t)  of  the  form 

n 

x { t )  =  I  a,  exp(X.  t) 
k=l  K  K 


and  details  of 
the  "complex 


(Bl) 


from  the  values  of  the  function.  The  function  x(t)  may  be  either  a 
scalar  or  vector  function  of  t  and,  accordingly,  the  ak  denote  either 
scalar  or  vector  constants.  We  consider  first  the  case  where  x(t)  and 
the  ak  are  scalars. 


Set 


nk  =  exp(Xkh) 

then  for  j  =  0,  ...,  n,  we  have  the  system  of  (n  +  1)  equations 

n  • 

x(t  +  jh)  =  akexp(Xkt) 

which  can  be  written  in  matrix  form  as 
x(t)  -  1  ...  -1 

x(t  +  h)  -  n1  ...  -nn 

x(t  +  nh)  -  n"  ...  -nn 


a1exp(A1t) 


anexp(Xnt)j 


(B2) 


(B3) 


(B4) 


Now  the  relation  expressed  by  Equation  B4  holds  if  and  only  if  the 
coefficient  matrix  is  singular,  that  is,  the  determinant  of  the  coefficient 
matrix  is  zero. 
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Let  A  +,  denote  the  cofactor  of  the  element  x(t  +jh).  Expanding  the 
0 

determinant  of  the  coefficient  matrix  in  terms  of  the  elements  of  the 
first  column  we  obtain. 

x(t+nh)  +  (An/Ap+1)  x(t+(n-l  )h)  +  ...  +  (A^A^)  x(t)  =  0  (B5) 

Hence,  a  function  x(t)  as  given  by  Equation  B1  satisfies  an  nth  order 
recursive  relation  or  difference  equation  of  the  form  given  by 
Equation  B5.  If  in  the  matrix  of  Equation  B4,  we  replace  x(t),  ..., 
x(t+nh)  by  1,  n...,  nn  respectively,  we  observe  that  the  coefficients 
^n+l-k^n+1  are  the  e1ementar-y  symmetric  (root)  functions  p,  (rh  .  ....  hn). 

pi  =  -(ni  +  •••  +v 


p2  =  nln2  +  ••*  +  ninn  +  n9n.,  +  •••  +  n,n„  +  ...  +  rj _ 


In  f,2M3 


n-1  n 


P„  =  (-1)  n,n2  ...  nn 


for  j=l  ,  . . . ,  n  and  set 


■  Wl 


xjk  =  x^j  +  ^ 


for  j-1 ,  ...,  n  and  k=l ,  ...,  n+1 .  Thus,  if  we  know  the  value  of  x ( t ) 
at  times  tj+(k-l)h,  then  from  Equation  B5  we  obtain  the  system  of  equations 


xn  •••  Xln  1  f  ci 


xnl  Xnmji  Cn 


for  the  coefficients  c. ,  ...,  c  . 

1  n 


ll  n+1 


n  n+1 


Suppose  Equation  B9  has  been  solved  for  the  coefficients  c^ ,  ...,  cR. 
We  can  determine  next  the  complex  frequencies  X^.  First,  one  solves  the 
polynomial  equation 


nn  +  c  nn"  1  +  ...,+  c.|  =  0 


(BIO) 


V 
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for  the  roots  n1 ,  • ...  Set  pk  =  |nkl  and  0k  =  arg  (nk)  where 
0  £  6k  <  2tr  then 

nk  *  exp[log  pk  +  1  (6fc  ±  2j7r)]  (811 ) 

and  Xk  is  given  by  the  equation 

hXk  =  log  pk  +  i(6k  ±  2j:r)  (B1 2) 

It  is  clear  from  Equation  B12  that  Ak  is  not  uniquely  determined  from  nk> 
However,  in  the  practical  situation  the  appropriate  value  of  Ak  is 
usually  clear. 

The  only  task  remaining  is  the  determination  of  the  amplitudes  ak. 
From  the  first  n  equations  of  the  system  given  by  Equation  B3  we  have 


’a-|exp{X^t)" 

a2exp(X2t) 


x(t+h) 


anexp^Xnt)J  |x(t+(n-l)h)J 


(B13) 


It  is  clear  that  this  linear  system  of  equations  is  solvable  for  the 
amplitudes  ak,  provided  the  coefficient  matrix  is  nonsingular.  However, 
the  particular  form  of  the  coefficient  matrix  lets  one  use  a  special 
method  for  solving  this  linear  system  of  equations.  Let  A^k  denote  the 
cofactor  of  the  (j,k)th  element  of  the  Vandermonde  determinant.  By 
Cramer's  Rule 


ajexp(Xjt) 


(  l  x(t+(k-l)h)A.  .)/(  I  n5"V .) 
k=l  kJ  k=l  J  kJ 


(B14) 


Dividing  both  numerator  and  denominator  of  the  right  hand  side  of  this 

equation  by  A  .  one  observes  that  the  ratios  A../A  .  are  the  elementary 
nj  kj  nj 

symmetric  (root)  functions  of  the  Hi/s,  excluding  n . .  This  observation 

K  J 

is  the  basis  of  the  algorithm  given  in  Reference  14  (page  275)  for 

computing  the  amplitudes  a..  In  order  to  see  clearly  how  this  algorithm 

J 

was  obtained  we  need  to  make  several  more  observations. 
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I 


Set 


then  the  denominator 


f (n)  =  (n-n-|)  •••  (n-nn) 
=r>n+Plnn_1+  •••  +  Pn 


(B15) 


n 


E  n 

k=l 


-  £  In-n, 


,  j  "kj  kn  dn  '  '  J 


f/(n-nj)  =  [f(n)-f(nj)]/(n-nj)  = 

(rin“ 1  +nn_ 2n .+  ...  *  m"-2+nf ’)  *  p,<n"‘W"V 


n-2,  n-3 


,  /  n-3^  n-4 

+Po(n  +n  n. 


‘ ‘  +nj"3^  +  ’ ' '  +p 


n-1 


(B16) 


+  n"'2)  (B17) 


Collecting  terms  in  like  powers  of  n  we  obtain  from  Equation  B17 
f/(n-nj)  =  nn  1  +  (r'j+Pi )in  2+(n2+p1nJ+P2)nn  3+(nj+P1n2+P2nj+P3)nn  4 
...  +(n'J"1+P1nj"2+p2nJ'3  +  ...  +Pn_-j ) 


+ 

(B18) 


k-1 

The  coefficient  of  n  in  Equation  B18  is  equal  to  the  ratio  A.  ./A  .  for 

n-k 

k=l  ,  ...,  n.  Hence,  if,  in  Equation  B18,  we  replace  n  by  x(t+(n-k)h) 

for  k=l ,  ...,  n  and  collect  like  terms  in  powers  of  n j >  then  we  can  express 

a.exp(A.t),  Equation  B14,  as  the  quotient  of  two  polynomials  evaluated 
J  J 

at  n j .  We  have  just  described  the  numerator  of  this  quotient.  The 
denominator  is  given  by  the  derivative  df/dn,  evaluated  at  n j .  This 
completes  the  derivation  of  the  algorithm  given  in  Reference  14  (page  275). 


It  is  clear  that  if  x(t)  is  a  vector  function  and  the  ak  are  vector 
constants  to  be  determined  we  could  apply  the  algorithm  described  above  to 
determine  the  jth  component  of  the  a^'s,  Tor  j=l ,  ...,  m.  That  is,  an 
equation  of  the  form  of  Equation  B13  has  to  be  solved  for  each  component. 
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Recently,  (References  2-4)  another  method  has  been  proposed  for 
determining  the  parameters  occurring  in  Equation  B1 .  Let  us  write 
Equation  813  as 


’x(t) 

-a,  ...  a„  - 

“exp(A.|t)“ 

x(t+h) 

= 

alnl  annn 

_x(t+(n-l)h) 

-  ainl  1  •”  annn  ■ 

.  exp(Xnt)_ 

(B19) 


Denote  the  column  vector  on  the  left  hand  side  of  Equation  B19  by  y(t), 
the  column  vector  on  the  right  hand  side  by  e(t)  and  the  matrix  by  A. 

The  Equation  B19  can  now  be  written  as 

y(t)  =  Ae(t)  (B20) 

Let  N  denote  the  diagonal  matrix  whose  diagonal  elements  are 
,  ....  nn-  Then 


and 


e(t+h)  *  Ne(t) 

y(t+h)  *  ANe(t) 


(B21 ) 

(B22) 


If  the  matrix  A  is  nonsingular.  Equation  B20  can  be  solved  for  e(t)  and 
we  have 


y(t+h)  =  ANA’^y(t) 


(B23) 


One  observes  the  Vandermonde  determinant  in  the  determinant  of  the  matrix 
in  Equation  B19.  Hence,  it  is  clear  that  A  is  nonsingular  if  the 
n-j ,  ...i  n  all  differ  from  one  another  and  zero  and  a^O  for  all  values 
of  k.  Note  that  Equation  B5  and  Equation  B23  are  equivalent  statements 
of  the  same  problem. 


Suppose  for  n  values  of  t,  t,  <  ...  <  t  the  vectors 

i  n 

yj  -  y(V 


(B24) 


V 
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and 


2j  =  y(tj+h) 


(B25) 


are  known.  From  Equations  B24  and  B20  we  obtain  the  matrix  equation 

[y1  •••  yn3  =  A[e1  •••  en]  ( B26 ) 

This  equation  can  be  written  symbolically  as 


$  *  AE 


(B27) 


In  a  similar  fashion  we  obtain  from  Equations  B25  and  B22 

y  =  ANE 

From  these  two  equations  we  obtain 

¥  $_1A  =  AN 

or,  on  setting  B  =  4>~^A 


H'B  =  4>BN 


(B28) 

(B29) 


(B30) 


Let  (uk)  denote  the  kth  column  of  A ( B )  then  we  can  replace 
Equations  B29  and  B30  by 

y*"1uk  =  Vk  (B31) 

and 

T  uR  =  nk$uk  (B32> 


for  k=l  ,  ...,n.  Thus  in  the  first  method  of  this  appendix  the  determi¬ 
nation  of  the  rt, ' s  was  formulated  as  a  problem  in  determining  the  roots 

of  a  polynomial  equation.  Equation  BIO.  In  this  second  method,  the  ru's 

-1  K 

are  displayed  as  eigenvalues  of  the  matrix  T<J>  or  of  the  generalized 
eigenvalue  problem.  Equation  B32.  Of  course,  here  too  the  determination 
of  the  nk's  could  be  reduced  to  the  determination  of  the  roots  of  a 
polynomial  equation. 


Examination  shows  that  the  two  methods  use,  or  can  use,  exactly  the 
same  data.  Since  the  eigenvectors  are  not  unique  we  have  no  assurance 
that  the  first  component  of  the  vector  satisfying  Equation  B31  is  the 
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ak  of  Equation  81.  It  frequently  suffices  to  know  the  ak  to  within  a 
multiplicative  constant.  Of  course,  once  the  nk's  are  known,  the  matrix 
E  is  known.  Hence,  if  E  is  invertible,  we  have  from  Equation  B26  an 
equation  equivalent  to  Equation  B13  and  thus  the  a^'s  are  determined  also. 


In  the  practical  situation  in  which  one  encounters  Equation  B1 ,  x(t) 
and  a^  are  m  dimensional  vectors  and  n  =  2m.  The  Equation  B1  can  be 
written  as 


x(t)  «  [a1  ...  an]e(t) 


(B33) 


Here  [a^  ...  a  ],  denotes  the  mxn  matrix  whose  columns  are  the  m 

dimensional  vectors  a^,  and  e(t)  is  the  n  dimensional  column  vector 

whose  jth  component  is  exp(X.t). 

J 


Again,  setting  y(t)  =  x(t+h)  we  have 

y(t)  =  [a1  ...  an]Ne(t) 


or 


y(t)  =  [n1a1  ...  nnan]e(t) 
Lastly  setting  z(t)  =  y(t+h),  we  have 

z(t)  =  [n]a1  ...  nnan]Ne(t) 

In  block  form  we  have  from  Equations  B33  and  B35 


rx(t,i .  r a'  -  a-i  e(t) 

[_y(t)J  j_n,a,  ...  VnJ 


and  from  Equations  B34  and  B35 


]  =  I"  3,1  1  Ne(t) 

J  IVl  '  * r*nan  J 


x(t+h)*|  ["  a1 

LyU+h) 


(B34) 


(B35) 


(B36) 


(B37) 


(B38) 


V 
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Denote  the  matrix  in  Equations  B37  and  B38  by  A.  If  A  is  nonsingular 
then  Equation  B37  can  be  solved  for  e(t)  and  we  have 


x(t+h) 


y(t+h) 


=  ANA 


(B39) 


Thus,  the  n-dimensional  vector  function 


z(t)  = 


(B40) 


satisfies  the  system  of  difference  relations  given  by  Equation  B39. 
Equation  B39  is  of  the  same  form  as  Equation  B23.  Hence,  as  we  have 
noted,  the  functions  z(t)  satisfying  Equation  B39  are  of  the  form 
uexp(At)  where  n  =  exp(Ah)  satisfies  the  condition 


and  the  vector  u  satisfies 


That  is. 


det[ANA  -nl]  =  0 


[ANA’  -nl]u  =  0 


ANA  =  nu 


(B41 ) 


(B42) 


(B43) 


n  is  an  eigenvalue  and  u  is  the  corresponding  eigenvector  of  the  matrix 
ANA'1  . 


The  matrix  ANA’  of  the  system  Equation  B39  is  not  known.  If  the 
vectors  x(t),  x(t+h)  and  x(t+2h)  are  known  for  n  values  of  t, 
t-j  <  ...  <  t  ,  then,  as  above,  we  can  determine  the  system  matrix 

-  _  1  -i  n 

ANA  =  W1  and  it  remains  to  solve  the  eigenvalue,  eigenvector  problem, 
either  Equation  B 31  or  Equation  B32.  It  is  clear  that  one  immediately 
obtains  the  eigenvectors  satisfying  Equation  B31  from  those  satisfying 
Equation  B32.  It  is  also  clear  that  the  eigenvectors  so  obtained  are 
not  necessarily  identical  to  the  columns  of  the  matrix  A  in  Equation  B37. 
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We  have  described  a  procedure  leading  to  an  eigenvalue-eigenvector 
problem  for  the  determination  of  the  complex  frequencies  when  x(t)  is  a 
scalar  and  also  when  x(t)  is  an  m-vector  and  n  =  2m.  This  same  procedure 
can  be  used  when  x(t)  is  a  2-dimensional  vector,  n  is  even  and  n  ^  1. 

The  details  are  so  similar  to  those  given  above  that  we  shall  not  repeat 
them  here. 

Examination  of  the  two  methods  of  this  appendix  shows  that  they  use, 
or  can  use,  exactly  the  same  data  in  the  scalar  case.  When  x(t)  is  an 
m-dimensional  vector,  the  first  method  requires  2n  readings  plus  n 

2 

reading  for  each  additional  component.  The  second  method  required  n 
readings  to  determine  4>  and  n  additional  readings  to  determine  V .  Thus 
when  x(t)  is  an  m-vector  the  two  methods  use  essentially  the  same  data 
also.  In  a  similar  fashion  we  note  that  it  takes  3n-2  readings  to  form 
4>  and  2  additional  readings  to  form  V  when  x(t)  is  a  2-dimensional  vector. 
Exactly  the  same  amount  of  data  is  required  by  the  first  method. 

If,  in  Equation  Bl ,  some  of  the  coefficients  ak  are  taken  to  be  zero, 
then  it  is  clear  that  the  determinant  of  the  matrix  4>  vanishes.  Or,  in 
other  words,  there  are  fewer  a^  and  to  be  determined  than  we  have 
allowed  for.  Thus,  the  vanishing  of  the  determinant  of  $  is  a  test  for 
determining  the  number  of  terms  in  Equation  Bl .  Evaluating  the  determinant 
of  a  matrix  is  usually  avoided,  if  possible,  and  in  numerical  work  a 
computed  zero  is  rather  rare.  Hence,  the  vanishing  of  the  determinant 
of  is  usually  not  a  good  test  for  determining  the  functional  form  of 
x(t) ,  Equation  Bl . 

The  process  of  determining  the  function  x(t).  Equation  Bl ,  is 
analogous  to  interpolation  with  trigonometric  polynomials,  in  the  sense 
that  the  coefficients  are  falsified  by  the  higher  harmonics.  The  number 
of  terms  in  Equation  Bl  is  not  specifically  known.  The  situation  may  be 
thought  of  as  follows.  Suppose 

x(t)  =  Z  a.exp(X.t)  +  z(t)  <844) 

k=l  K 
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There  is  nothing  which  prevents  one  from  forming  the  matrices  $  and  ¥ 

(or  determining  the  coefficients  c^  ....  cn,  Equation  B9,  as  described 
above.  However,  it  is  clear  that  these  quantities  will  be  falsified  by 
the  contribution  of  z(t).  Equation  B44  also  indicates  why  the  vanishing 
of  the  determinant  of  4>  is  not  a  clear  cut  process.  Taking  a  larger 
number  of  terms  in  Equation  B1  might  possibly  reduce  the  influence  of  z(t) 
on  the  parameters  A ^  and  a^  of  Equation  B1  . 

The  columns  of  the  matrix  A  defined  by  Equation  B1 9  for  the  scalar 
case  and  by  Equation  B37  for  the  case  where  the  a^  are  m  dimensional 
vectors,  represent  the  eigenvectors  which  are  to  be  determined.  Note 
that  these  eigenvectors  have  a  somewhat  restricted  form.  That  is,  for 
the  scalar  case,  the  second  entry  in  the  k  th  column  is  times  the 
first  entry.  Similarly,  for  the  m  dimensional  case  the  last  m  entries 
in  the  k  th  column  are  just  the  first  m  entries  multiplied  by  n^. 

This  special  form  of  the  eigenvectors  is  the  basis  for  the  modal 
confidence  factor  introduced  in  (Reference  5).  Thus,  how  well  the 
eigenvalue  and  eigenvector  solutions  of  Equation  B31  satisfy  the  special 
form  determines  the  degree  of  confidence  that  the  solutions  so  determined 
represent  the  complex  modes  and  frequencies  of  the  structure.  The  modal 
confidence  factor  can  be  used  in  addition  to  or  as  an  alternative  to  the 
ratio  of  successive  determinants  test. 

We  conclude  this  section  with  an  examination  of  the  matrices  A 
defined  by  Equation  B37  and  E  defined  by  Equations  B26  and  B27.  If 
tj  =  t^  +  (j-l)h  for  j  =  2,  ....  n  then  the  matrix  E  can  be  written  as 


_  ~ 

r  n-1  1 

E  = 

expA-|  t 

1  n,  ...  n, 

•  •  • 

•  •  • 

expApt^ 

n-1 

V  nn  nn  j 

Clearly  E  is  nonsingular  if  the  n-j ,  ....  np  are  all  different  from  one 

another.  Now  =  exp  (A^h).  Let  us  write  A^  =  +  if^,  and  suppose, 

for  some  k  and  j^k,  nt/ri.  =  exp  h(A,-A.)  =  1. 

*  J  K  J 
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This  last  equality  will  hold  if  and  )h(8^-8j)|  *  2ul.  Thus 

even  though  the  ,  ...  Am  are  all  different  there  is  the  possibility  of 
a  sampling  rate  for  which  the  matrix  E  is  singular.  It  seems  readily 
apparent  that  there  should  be  no  difficulty  in  choosing  a  sampling  rate 
which  will  insure  that  the  matrix  E  is  nonsingular. 

The  columns  of  the  matrix  A,  Equation  B37,  represent  the  eigen¬ 
vectors  associated  with  a  linear  system  of  first  order  difference 
equations.  Accordingly  the  matrix  A  should  be  nonsingular.  The 
function  x(t).  Equation  B33,  perhaps  more  precisely,  should  be  written 
as  x(t)  =  [c-j  a1  ...,  cnan]e(t).  In  Equation  B33  the  coefficients  c^  are 
included  in  the  symbol  a^.  Thus  the  function  x(t)  is  a  linear  combination 
of  the  eigenfunctions  a^expA^t.  It  is  clear  that  in  a  particular 
experiment  one  or  more  of  the  eigenfrequencies  A^  need  not  be  excited. 

That  is,  those  columns  of  A  corresponding  to  the  nonexcited  Afc  are 
columns  of  zeros. 

From  the  point  of  view  of  determining  the  matrix  A  and  the  eigen¬ 
frequencies  Ak  it  follows  from  our  remarks  that  the  set  of  sampled  values 
of  the  function  x(t)  may  not  be  adequate  for  determining  all  the 
eigenvectors  and  frequencies.  On  the  other  hand  there  is  the  possibility 
that  one  has  allowed  for  more  eigenfrequencies  in  some  selected  frequency 
range  than  are  physically  present. 

The  difficulties  noted  were  recognized  by  S.  R.  Ibrahim,  the  principal 
author  of  the  method  being  discussed,  and  were  treated  in  his  papers 
(References  3  and  4). 
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APPENDIX  C 

SOME  MATRIX  EIGENVALUE-EIGENVECTOR  CONSIDERATIONS 

In  this  section  we  present  and  examine  material  which  we  feel  is  the 
basis  for  the  method  of  Wittmeyer  (References  6  and  7).  We  begin  by 
considering  the  equation 

[-u)2  A  +  B]y  =  r  (Cl) 

A  and  B  denote  symmetric  matrices  of  order  m.  The  matrix  A  is  real 
while  B  is  generally  complex.  An  equation  of  this  form  represents  the 
dynamic  equation  for  a  system  with  m  degrees  of  freedom  and  structural 
damping. 
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We  take  as  our  first  task  the  determination  of  the  eigenvalues  A^ 

and  their  corresponding  eigenvectors  u^.  In  the  practical  situation 

the  known  data  are  the  values  assigned  to  w  and  the  excitation  vector  r, 

and  the  observed  responses  y.  However,  for  the  present,  we  suppose  the 

matrices  A  and  B  are  known  also.  The  starting  point  for  the  determination 

of  the  eigenvalues  and  eigenvectors  is  the  observation  that  for  an 

excitation  vector  r  =  r.  satisfying  the  condition 

<3 

Vj  =0  (C8) 

2 

for  all  k  i  j  and  for  any  ui  /  for  all  k,  the  solution  y  =  y.  is 

K  J 

*5  ’  uj(ujrj)/mj(xj  '  “2>  <C9) 

That  is,  y.  differs  from  the  eigenvector  u.  by  the  scalar  factor 

T  J  2  J 

u.r./m.U.  -  »  )• 

It  is  rather  unlikely  that  an  r  such  as  r.  is  known.  However,  if 

J  p 

A.  is  well  separated  from  the  A,  ,  for  k  ?  j  and  if  u  is  "close"  to  A., 

J  ^  3 

it  is  clear  from  Equation  C6  that  u.  is  the  principal  contributor  to  the 

2  J 

value  of  y.  Thus  for  w  close  to  A .  we  may  suppose  the  initial  guess 

J 

of  r  for  r j  expressed  in  the  form 


r  =  I  p  Au. 
k=l  k  k 


(CIO) 


and  in  particular  p.  f  0.  Then  by  Equation  C6,  the  corresponding  y  is 
J 


in  2 

y  =  Z^  P|<u|/^k  ”  ) 


(Cll) 


If  the  matrix  A  were  known,  we  could  iterate  on  the  two  operations 
represented  by  Equations  CIO  and  Cll.  That  is,  we  take 


r^  *  Ay  =  Z  P kAU|c/(X|C  -  w2) 


(Cl  2) 


y(1)  -  Pkuk/(Ak  -  O)2)2 


(Cl  3) 
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and  so  on  to  obtain 


»  - 


—  E  pkuk^k  ^  ^ 


2\n+l 


k=l 


(Cl  4) 


It  is  clear  from  Equation  Cl 4  that  if  for  each  value  of  n,  y^  is 
normalized  then 


y(n)  _  y  (n)  +  const,  u. 

J  J 

Thus  for  n  sufficiently  large,  we  will  have 


(Cl  5) 


=  Ay.<">  .  d-Au- 
J  .1  3  3 


(n) 


then  from  this  r.  we  obtain 


yt("+1)  ’  ^juj/(xj  •  »2) 

and  from  r^n+^  =  Ay^n+^  we  obtain 
3  3 

y(n;2>s  -  J)2 

We  have  immediately  from  the  ratio  of  the  vth  component  y^ 
vth  component  of  y/n+^;  that  is 


(Cl  6) 


(Cl  7) 


(n+1 ) 


to  the 


^Vty^VV*2  lC’8) 

and  so  A.  is  known  also.  This  is  one  way  in  which  the  eigenvectors  u. 

3  3 

and  the  corresponding  eigenvalues  A  ^  can  be  determined. 


We  are  going  to  describe  now  yet  another  way  for  determining  the 

eigenvector  u.  and  the  corresponding  A..  A  given  vector  y  may  be 
3  3 ' 

represented  as  a  linear  combination  of  the  eigenvectors  u^.  Thus 


m 

y  =  T  nku. 
k=l  K  K 


(Cl  9) 
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«  ' 

! 


Then  we  readily  establish  that 


r  =  [-</  A  +  B]  y 


m 


(C20) 


=  E  nk( Xk-ai  )  Ail,, 


k=l 


For  this  second  method  we  start  with  an  initial  excitation  vector 
*  and  two  frequencie 
represented  in  the  form 


r^  and  two  frequencies,  w2  and  to2  close  to  A..  Again  supposing  r^ 

J 


we  have 


=  Z  okAuk 


-  [-„2  A  ♦  B]-1  r<°> 

m  2 

~  E  pkuk//^k"w  ^ 


(C21) 


by  Equation  C6.  Take 


r<°>  *  [-u>2  A  *  B]  y(0) 

,  ~2 
Ak  -  u) 


m 

"*  E  P. 


.  ,  'k  .  T  auk 
k=l  A.  -  w 


(C22) 


by  Equation  C20.  Now  take 

.0)  .  _(o)  .  p(0) 


m  -2 

-  i  SL- 


■K  Pkauk 


(C23) 


k=l  A^  -  w 

Considering  the  computations  indicated  by  Equations  C21-C23  as  an 
iteration  step,  we  have  after  (n+1)  steps 


#  _  \  ro  / ~  2  2 ,  n 

y  '  1  ^  n+T  pkuk 


(C24) 


k*l  (Ak  -  u)c)' 
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If  the  A^  are  well  separated  and  if 

~2  2 
U)  -  0) 


A  .  -  co 

J 


=  1 


(C25) 


,(n) 


then  obviously  y  '  ->  const,  u.  and  A.  can  be  determined  in  the  same 

J  0 

fashion  as  above.  Equation  C18.  In  general,  the  equality  expressed  by 
Equation  C25  is  not  satisfied.  In  this  event  some  normalization  of  y^ 
will  probably  be  required  with  each  iteration  step.  The  results  obtained 
by  assuming  Equation  C25  held  are  not  altered  however. 


In  the  methods  discussed  in  this  appendix  up  to  this  point  we  have 

assumed  that  the  eigenvalues  are  well  separated.  Now  we  consider  the 

case  of  two  close  eigenvalues,  say  Aj  and  A .+.| .  We  consider,  again,  a 

constant  vector  r.  satisfying  the  condition  uJr .  =  0  for  all  k  f  j. 

J  K  J 


Thus  if  we  excite  with  r.  and  frequency  u>.  close  to  A.  we  obtain 

J  J  J 


(C26) 


Similarly,  if  we  excited  with  the  same  r.  and  a  frequency  w.  we  obtain 

*3  J 


y.  =  (UjVjlKj/mjUj  -  „.2) 


From  Equations  C26  and  C27  we  obtain 


yj  ‘  yj  (xj  •  “j2)/(xj  -  “/> =  yj  -  wj  * 0 


(C27) 


( C28) 


Now  suppose  we  have  linearly  independent  vectors  rg  and  r^  and  that 

r.  =  ar_  +  br. 

J  a  b  ( C29) 

It  follows,  because  of  Equation  C6,  that  the  excitation  vectors  and 

A  /v 

frequencies  r  ,  w.,  r,  ,  w.  and  r,  ,  u>.  will  determine  response  vectors 
a  ®  AJ  ^  J  Uj 

y  ,  y.,  y  ,  and  y.  respectively.  It  follows  also,  from  Equation  C29,  that 

3  D  3  D 


yj  =  aya  +  byb 


(C30) 


yj  =  aya  +  byb 
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From  Equation 

C28  we  have  the  relation 

*  cyj  ' 

tv,b3[;] 

+  P^a’  yb] 

= 

[».•  yt]  + 

p[ya>  ybl] 

from  which  we 

have  finally 

tV][v 

*  p!V 

£ya’  yb]  f 
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APPENDIX  D 

DETERMINATION  OF  MASS,  DAMPING  AND  STIFFNESS  MATRICES 

In  this  section  we  want  to  consider  the  problem  of  determining  the 
mass,  damping  and  stiffness  matrices  directly  from  experimental  data. 
However,  first,  for  completeness  and  understanding  we  examine  the 
problem  of  solving  systems  of  linear  equations,  the  idea  of  a  generalized 
solution  and  of  a  pseudo  or  generalized  inverse. 


Let  A  denote  a  real  symmetric  matrix  of  order  m.  For  each  such 

matrix  A  there  is  an  associated  set  of  m  real  numbers,  A,,  ...,  A  and 

1  m 

a  linearly  independent  set  of  vectors  u^ ,  u^  which  satisfy  the 
condi tions 


Auk '  Vk 


u]a  .  Xju] 


(Dl) 


and 


If  we  set 


U  =  [uj  ...  um]  (D2) 

then  clearly 

UTU  =  I  =•  UUT  (D3) 


and  the  matrix  A  can  be  written  in  the  form 

A  =  UAUT 


(D4) 


where  A  is  a  diagonal  matrix  with  the  A^’s  along  the  main  diagonal. 


The  properties  of  the  matrix  A  as  an  operator  are  readily  apparent 
when  A  is  expressed  in  dyadic  form.  Thus,  if  we  consider  Ax,  then  we 
obtain  from  Equation  D4  the  following  representation  for  the  matrix  A 


A  “  A-j  Ui  u-| 


m  m  m 


(D5) 
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From  F.quation  05  it  is  clear  that  the  range  of  A  is  the  space  spanned  by 
those  vectors  uk  for  which  Ak  t  0.  It  is  clear  also  that  A  is  1-1  on 
the  range  of  A.  Hence,  if  the  range  and  domain  of  A  are  the  same,  that 
is,  if  ^  0  for  all  k,  then  A~^  exists  and 

A'1  =  (1/A-j  )u-ju|  +  ...  +  0/Am)u1u{  (D6) 

Let  denote  the  set  of  values  of  k  for  which  A^  f  0  and  Ng  the 
set  of  values  of  k  for  which  A^  =  0.  From  our  remarks  above,  Ax  =  b 
has  a  "strict"  solution  if  and  only  if  b  lies  in  the  range  of  A.  Now 
b  lies  in  the  range  of  A  if  and  only  if 

ujjb  =  0  (07) 


for  all  k  in  Ng.  Thus  if  b  satisfies  these  solvability  conditions  then 
a  solution  to  Ax  =  b  is  given  by 


x  =  l  (1/Aju.iJb 
keN-j 


(08) 


If  we  set 


y  =  x  + 


Vk 


(D9) 


then  Ay  =  b  also  and 

I  |y|  I  >  IN 


(010) 


Even  if  b  does  not  satisfy  the  solvability  conditions,  we  can 
still  compute  an  x  by  means  of  Equation  D8.  Then 

Ax  =  l  u.u^b  (Dll) 

keN1  K  K 

and 


b  -  Ax 


(012) 
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Then  for  any  y  whatsoever 

b  -  Ay  =  E  c.u,  +  E  u.ulb  ( D1 3) 

keN-j  K  K  keNQ  K  K 

It  follows  from  Equation  D12  and  D13  that 

I  |  b  -  Ay  ||  _>  |  |  b  -  Ax||  (D14) 


Set 


a1  =  e  d/uv; 

keN}  K  K 


(D15) 


A*  is  called  the  generalized  or  pseudo  inverse  of  A.  We  have  seen  that 

if  X.  is  not  zero  for  any  value  of  k  then  A^  =  A  ^ .  If  the  equation 
k  I 

Ax  =  b  has  solutions  then  x  =  A  b  is  the  solution  of  least  magnitude. 
Equation  DIO.  On  the  other  hand,  if  Ax  =  b  has  no  solutions  then  x  =  A*b 
satisfies  the  equation  Ax  =  b  as  well  as  or  better  than  any  other  vector  y 
(Equation  D14). 


We  have  characteri zed  the  generalized  inverse  for  a  real  symmetric 
matrix  A.  Next  we  want  to  characterize  the  generalized  inverse  of  an 
arbitrary  real  (m  x  n)  matrix  A,  where  01  ^  n,  and  of  its  transpose  A^. 

To  achieve  our  purpose  we  consider  the  matrix  A^A. 

Now  the  matrix  A^A  is  real  symmetric  and  of  order  n.  Hence  there 
are  real  eigenvalues  ,  ...,  f  and  a  complete  orthonormal  set  of  n 

dimensional  eigenvectors  u-j  ,  ...,  u  .  It  follows  from  the  equation 

■  Vk  ’  5k  (tll6) 


that  2  Set  xk  =  1  0  and  set 


vk  a  OAk)Auk,  if  Xk  t  0 


vk =  Vk =  0 


(D17) 


if  \k  =  0.  One  readily  shows  that  the  vk,  for  k  in  Np  constitute  a 
set  of  orthonormal  m  -  dimensional  vectors. 
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We  obtain  from  Equation  D17,  regardless  of  the  value  of  A^, 
Atu^...,  U^]  -  [vi  »  •••»  V^]A 

From  this  equation  we  infer  that  A  has  a  representation 

A  =  Z  A.  v  .  uT 
keN7  K  K 


(D18) 


(D19) 


then,  as  above,  we  take 


keN 


We  have  from  Equation  D1 9  that 


A  M  <1/xk>Vk 


AT  =  Z  A.u.vJ 
keN1  K  K  K 


and  hence 


T  1  T 

(A1)  =  Z  (1/A  )  V.  u/ 

keN7  K  K  K 

In  a  manner  similar  to  the  above  one  can  obtain  results  of  the  same 
nature  for  complex  matrices. 


(020) 


(021) 


(D22) 


We  now  describe  a  rather  simple  and  obvious  approach  for  determining 
the  mass,  damping  and  stiffness  matrices.  In  general  this  procedure 
fails  in  practical  situations  and  hence  is  really  no  method. 

The  steady  state  response  of  the  system  of  equations  Mx  +  Cx  +  Kx  =  f 
to  a  harmonic  excitation  f(t)  =  r  exp  (iu^t),  r  a  constant  vector,  is  of 
the  form  x(t)  =  y  exp(itu.|t)  where  y  is  a  constant  vector.  If  for  k  =  1, 
...,  m  we  take  r  equal  to  the  unit  Cartesian  coordinate  vectors 
respectively  we  obtain  the  matrix  equation 

M  +  K  +  iu)jC]  [y^ ,  ...»  ym]  =  I  (D23) 


~0 


V 


A 


mm  ~  -  a- a 
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Since  u>i  is  not  a  characteristic  value  we  can  infer  that  the  matrix 
[y^ ,  ....  ym]  is  nonsingular.  It  follows  then  from  Equation  D23  that 


M  +  K  =  Re[y, ,  ....  Ym] 


j-jC  =  Im[y-j ,  •••»  yml 


(D24) 


(D25) 


If  we  repeated  the  process  for  w  =  then  we  would  have  another  equation 
of  the  form  of  Equation  024,  say 


-u)2  M  +  K  =  Re[z-j . zm]_1 


( D26) 


The  Equation  D25  gives  the  damping  matrix  C  and  Equations  D24  and  D26 

determine  the  mass  and  stiffness  matrices.  It  comes  as  somewhat  of  a 

surprise  that  the  M,  C  and  K  determined  by  this  method  are  not  very  good. 

Of  course,  one  immediately  concludes  that  the  matrices  [y-j ,  ...,  y  ]  and 

[z, ,  ...,  z  ]  must  be  ill  conditioned.  If  one  looks  at  Equation  A34  or 

Equation  A38  of  Appendix  A  the  reason  for  the  ill  conditioning  becomes 

apparent.  Thus  /<>-|  will  be  close  to  X.  say  and  for  k  =  1,  ...,  m  the 

y,  will  have  nearly  the  same  direction,  as  the  modal  vector  u..  This 
K  J 

observation  coupled  with  the  limited  accuracy  of  experimental  data 
explains  the  poor  results  obtained  by  this  method. 

It  is  clear  that  one  can  obtain  a  better  set  of  spanning  vectors 
[y-j  >  . y  1  by  a  judicious  choice  of  frequencies,  ,  ...,  wm. 

However,  the  resulting  set  of  equations  are  more  complicated.  This 
method  was  discussed  recently  in  Reference  8  and  we  want  to  examine  the 
method  here. 


In  Reference  8  it  was  assumed  that  the  mass  matrix  M  was  known.  It 
is  instructive  to  consider  first  that  all  three  matrices  M,  C  and  K  are 
unknown. 
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Let  yk  exp  (iwkt)  denote  the  steady  state  response  to  the  harmonic 
excitation  rk  exp  (iu^t),  rk  real.  We  have  then  the  equation 

-"“kJ'k  *  ^k +  ,cVk  ■  rk  (°27) 

Set  yk  =  yRk  +  iy^.  Then  Equation  D27  can  be  rewritten  as  two  real 
equations 


+  KyRk  +  C("^yT^  =  r 


kJIk'  ’k 


M('“kyIk)  +  KyIk  +  C(ukyRk^  =  0 


(D28) 


Now  let  us  suppose  m,  the  order  of  the  matrices  M,  C  and  K  is  even 
and  set  p  =  3m/2.  Also,  set 


~T  _  r  2  T  T  T  , 

yk  “  ^VRk’  yRk*  "Vlk^ 

yp+k  =  ^wkyIk’  yIk’  wkyRk^ 


(029) 


For  k=l  ,  p  we  have  then  the  system  of  equations 

[M  K  C][y-| ,  ....  yp  yp+1,  ....  y2p]  =  [r^ . rp  0  ...  0]  (030) 

A  A 

where  [M  K  C]  is  an  m  x  3  m  matrix,  [y  ...  y2p]  is  3m  x  3m  and  [r^  ...  rp, 
0  ...  0]  is  m  x  3m.  In  view  of  our  remarks  above  on  generalized  inverses 
we  may  write 

[M  K  C]  =  [r1  ...  rp  0  ...  0]^  ...  y2p] 

where  if  [y1  ...  y2p]  is  not  strictly  invertible,  then  [y^  ...  y2p]-1 
denotes  the  appropriate  generalized  inverse  of  [y-j  ...  y2  ]. 


Let  us  observe  that  if  any  row  of  the  matrix  [r^  . . .  rp  0  ...0] 
consists  entirely  of  zeros  then  the  same  row  of  the  matrix  [M  K  C]  will 
have  only  zeros.  It  follows  that  the  r^ ,  ....  rp  should  span  an  m 
dimensional  space.  Alternatively,  it  follows  that  the  system  determined 
by  the  matrices  M,  C  and  K  should  be  excited  at  least  once  at  each  station, 
if  we  are  to  avoid  a  row  of  zeros  in  the  matrices  M,  C  and  K. 
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For  the  case  where  the  matrices  M,  C  and  K  are  of  order  m  (i.e.  for 
a  true  m  degree  of  freedom  system)  it  is  possible  to  choose  the  amplitudes 
r^  and  the  frequencies  so  that  the  matrix  [y^  ,  ....  y2^]  has  a  strict 
inverse.  That  is,  in  principle  it  is  possible  to  recover  the  matrices  M, 

C  and  K  from  experiments.  For  the  practical  case,  however,  many  of  the 
points  at  which  one  observes  the  response  to  an  excitation  are  not 
suitable  points  for  exciting  the  structure. 


For  the  reasons  just  mentioned  determining  the  matrices  M,  C  and  K 
by  the  method  represented  by  Equation  D31  has  very  limited  applicability. 
The  situation  is  quite  different  if  the  mass  matrix  M  is  known.  In  this 
case  set 


and 


-T  r  T  T  n 

yk  "  tyRk’  'Vlk-* 

ym+k  =  ty{k’  wkyRk^ 


(D32) 


?k  ‘  rk +  "k%k 
rm+k  '  “kMyIk 


Then  for  k=l ,  . . . ,  m  we  have  the  system  of  equations  for  the  m  x  2m 
matrix  [K  C] 

Ck  c3 [y-j »  ...»  y2m]  =  Cr-]>  r2m]  (D33) 

2 

The  vector  r^  in  the  expression  r^  +  w^My^  can  be  the  same  vector  for 
all  values  of  k.  That  is,  if  the  matrix  M  is  known  then  exciting  the 
system  at  a  single  well  chosen  station  with  a  set  of  well  chosen 
frequencies  ,  . . . ,  the  matrix  [y^ ,  ...,  y^]  will  be  strictly 
invertible  when  the  eigenvalues  A^  are  well  separated. 

If  for  some  value  of  k,  A,  is  an  eigenvalue  of  multiplicity  greater 

/\  A  _ 

than  1  then  it  is  clear  from  Equation  A38  that  the  matrix  [y^ ,  .  ..y^J, 
obtained  by  exciting  the  system  at  a  single  point,  will  be  singular. 
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Suppose  the  eigenvalue  A^  is  of  multiplicity  p  and  for  simplicity 

suppose  also  that  A^  is  the  only  multiple  eigenvalue.  In  this  case  we 

need  p  suitable  points  at  which  to  excite  the  system.  (Even  though  p 

excitation  points  are  required,  the  system  is  excited  at  only  one  point 

at  a  time.)  Then  for  near  resonance  frequencies  a>,  ,  .  .  .  ,  co  ,,  we 

i  m-p+i 

should,  at  least  in  principle,  obtain  a  matrix  [y^ ,  ....  y2m]  which  is 
strictly  invertible.  For  the  frequency  which  is  close  to  the 
multiple  eigenvalue  A^  the  system  should  be  excited  at  each  of  the 
p  excitation  points,  one  point  at  a  time. 

Thus  corresponding  to  each  multiple  eigenvalue  we  need  a  set  of 
excitation  points  equal  in  number  to  the  multiplicity  of  the  eigenvalue. 
The  same  excitation  points,  if  suitable,  may  be  used  for  different 
eigenvalues.  We  believe  a  similar  procedure  shoud  be  followed  for  close 
eigenvalues  also. 


Set  YR  =  ^Rl’  **•  YRM-'1  YI  ^11  ’  YIM^’  R  '  £r1*  rm-'  and 

let  Q  denote  the  diagonal  matrix  with  diagonal  elements  for  k=l ,  ...,  m. 

Then  Equation  D33  can  be  written  in  block  matrix  form  as 


KYr  -  CYjfi  =  R  +  MYRfi2 
KYj  +  CYRft  =  MYjfi2 


(034) 


If  we  multiply  Equations  D34  (on  the  left)  by  M  ^  we  obtain  a  system  of 
equations  for  the  transformed  stiffness  matrix  K*  =  M"  K  and  damping 
matrix  C*  =  M~^C,  namely 


K*Yr  -  C*Yjfi  =  R*  +  YRn2  (035) 

K*Yj  +  *CYRn  =  Yjf22 
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If,  for  example,  we  multiply  the  second  of  the  Equations  D35  on  the 
right  by  Yj_1Yr  and  subtract  the  first  equation  therefrom  we  obtain 

C*CYRnYI“1 Yr  +  Yjft]  =  Yjft2  -  R*  -  YRn2  (p36) 

Once  this  equation  is  solved  for  C*  the  value  so  obtained  can  be  used 
in  either  of  the  Equations  D35  to  determine  K*.  A  procedure  of  this  kind 
is  probably  preferable  for  determining  the  matrices  K  and  C  or  their 
transforms  K*  and  C*,  to  solving  Equation  D33  for  the  matrix  [K  C],  for 
example. 

We  have  examined  the  method  (Reference  8)  primarily  from  the  point  of 
view  of  an  evenly  determined  system.  If  the  systems  Equations  D33, 
Equation  D34  or  Equations  035  were  over  or  under  determined  then  the 
methods  described  at  the  beginning  of  this  appendix  could  be  used  to 
determine  matrices  K(K*)  and  C( C* ) .  It  should  be  clear  from  our  remarks 
above  that  even  though  the  number  of  experiments  exceeds  the  number  of 
unknowns  we  do  not  have,  necessarily,  an  over  determined  system.  Once 
the  matrices  K(K*)  and  C(C*}  are  determined  one  can  calculate  other 
quantities  of  interest  as  described  in  Reference  8. 
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Figure  1.  An  Example  of  an  Element  of  the  Frequency 
Response  Function  Matrix 

Graph  in  a  Comdex  plane  of  the  function 

6  _ 

w(u)  =  e(w)  +  in(w)  =  I  [bk/(iw-Xk)  +  ( ico-Xk ) 

k  “  1 

The  characteristic  values  X^  and  the  coefficients  b^  are  piven  in 
Table  1A.  The  boxes  are  the  values  of  w (to)  for  the  frequencies  listed 
in  Table  1A. 
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B 


12 


f  the  Argument  of  w(u>) 


=  5.2 

u  =  4.7 

x  Plane  of  £(u>)  +  in(w)  =  w(w)  -  b^/ 
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10 


B 


H 


! 


2 


H - 1 - 1 - 1 - (-*•  U) 

H  4.7  £  5.2  E 


Figure  5. 


"Upper  Curve"  -  Graph  of  |w(co)|  "Lower  Curve"  - 
Graph  of  |w(w)  -  t>4/(iu)  -  X^)| 


V 
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b4  =  0.866025404  -  0.5i 


Graphs  of  w(w)  for  Four  Values  of  b. 
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TABLE  I A 


DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED 
t>4  =  0.866025404  +  0.5i 


CHflRflCTERISTIC  VflLUES 

Ak 

LTeflflflflQQflP  fti  4  *vw vwwVW.rtft  _Q  eaQQAOQaArjn 

r  JwKmXKIC  iu  X.  O.  '^wwww.  uc 

tneftftVMww  fti  4  fywvyyyyvy « 4  MQQAQQQfir. Ai 

djWWwix'in.  i.  uuwuwwxtu 

7  /MMMWwvMVir  A4  O  flQflQQQflftflr  fti 

"i.  WwWWW3L"gl  7.  MMWUAL'w  c  wvwwOK  vi 


COEFFICIENTS  bk 

-L  736481780E-01 
9. 396926210E-01 
-7.  660444430E-01 

FREQUENCIES  u 

t/vvwvvvvvv- 1 

WKwAaWA.  ’  w 

2^WWWWW  .  ftp 

3<yvwvvvwyy ,  m 

4  Tftflflftflflftflr  .IMI 
T.  I  WVWW^MW 

8rwwwww . 

jwwwJw^tw 

L  oJWWUWXTTli 

FUNCTION  VALUES 

-4.  292738520E+00 
-3  24913691 5E+08 
4.  389415898E+00 
-2  339350993E-01 
-1 669958191E+00 
-4  715254288E+00 
5. 997219248E-01 
7. 614061425E+00 
3. 7774G4651E*# 

4. 038931343E+00 
-2  4404U390E+60 
-1  714513046E+00 


9  848077538E-01 
-3.  420201430E-01 
-6. 427876100E-01 

oj, 

6+p 

4  ■XMVMWWWir.flfl 

i.  jWWWWL'w 

3i#W WVWWVW.ftfl 

.  lWWWWCTw 

4  4  fMVMMVWW.M 

5'VWWVWW ,  nr* 

■  tWWWAAA.  T  w 
.  pWOWwwltw 

l  iToeeeeeeE+ei 


2  251750778E+00 
7. 6277632 79E+00 
2  751933130E+00 
1  165997917E+00 
-3.  623276859E+00 
-5.  223206098E-01 
5. 073768837E+00 
3.  842118338E+00 
-4.  898518621E+00 
-3. 820507362E+00 
-7.  779925282E-01 
1  169627699E+00 


3. 420201430E-01 
8.  660254040E-01 
-9  396926210E-01 


TRUNCATED 

4  |VUVVMVMMriM 
7  n44MMMMWI4ir.«4t 

A.  iWWWWX'w 

.  jwuwww>'m,! 

7  vwwwyyyy  fu 
6  iwwwuw.^1 

fOQflQflflflCjJfl 

WWWWWV’W 

-4.  710000000E+00 

5/VWVWUVMW  fti 

,  ywwwwuTIJ, 

3.  770000000E+00 
4.  030000000E*00 

_7  4  4IMWMWIflrtT.flfl 

C.  WwwwwXTw 

-1  710000000E400 


VALUES 


2AAQQAQAQAr  >M 

5.  000000000E 1 00 

L4  ft/WWWVWWir  .  AJ 


9  396926210E-01 

SMWWtftfMWMP  Ai 

.  wl 

1 420O1438M1 


VALUES 

C. 

i .  Ocwwvw.  MW 

C.  I  JVWWVVW  ■  IW 

t4  /TMMWMMlf  .Art 

IwwwwwX’w 

JfVMMMlflflflr.flfl 
.  Dtuwww».Tw 

*5. 

5ftTft<WWIft<W  .  ftft 

3n4fWMMMMW  .  <M> 
4  AQAAQAQAQPiM 

“  OJWWwW.’W 

7  AV^MAMTiM 

y 

7'VWVWVWVW  fti 
,  #1, 
7  4  flvvwwr.^ 

y  XD 
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TABLE  IB 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION 


NUWER  OF  ITERATIONS*  12 


RADIUS  OF  DIFFERENCES  3.  6411826488-08 


COfFUTED  CHARACTERISTIC  VALUES 

-1  7499999908-01 
-12590008408-01 

-2  9999999798-01 


4  OAaQOQQIKriiU 

1. 

3.  9999999988+00 

o  QottkNkmr  ■  />ft 
y.  iwwwvfLfW 


ft  CQflQQflMW  ftft 
O.  ’  Ow 

-1  0000008388-01 

_ft  QMAftOAi  nC—A4 

L. 


CHARACTERISTIC  VALUE  DIFFERENCES 
-L  0350600008-09  -5. 

4  003000000E-09  1 

— O  /vvvvvvvvv‘  <vi  -j 

c.  tJWW3WC"t!7  “i. 


ftftrtflflflftftftr  Aft 

/-liVWMWWV-  ftft 


ft  ftMUVVVVW^  i  i 
7  ftftftft/MWw  ftft 

3.  W— R? 

1  226000690E-09 


COMPUTED  COEFFICIENTS 

-1  736481367E-01 
9. 396926567E-01 
-7.  66M44440E-01 


9  8480775398-01 
-3.  4282013888-01 
-6.  4278758938-01 


3. 428201709E-01 
8.  669254162E-01 
-9  3969262698-01 


COEFFICIENT  DIFFERENCES 


■4. 129180000E-88 
-3  5689080088-08 
9  770000008E-10 


-8  5800690088-10 
-4  219000600E-09 
-2  071700800E-08 


-2  702200000E-08 
-1  224000000E-08 

5iv/whww  /V> 


VALUES 


.  UUWWUUx^w 
.  777777^+W 

11000000608+01 


-2  6300000008-09 

4iftivwvwvw  aft 
i  mmmaaat  ift 


9  3969261358-91 

SflflOftQftiftflf  fli 
.  WWWlwCDl 

3  420201491E-01 


7  5240600008-09 
-18921000008-08 

ftrtftftftftftftnr  flft 
v. 
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TABLE  1C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS^  12 

RADIUS  OF  DIFFERENCES  1  62O600949E-fte 

COMPUTED  CHARACTERISTIC  VflLUES 

X, 

-i  752927B73S-m  k  1 200003610E+B0  -ft  507245B55E-K  10M099M1E«0 

-1  248179243E-01  3. 9999ft0107E+€0  -9. 9ftl410287E-€2  5. 000273706E*« 

-2. 9687Uft76E-81  9. 001500881E+60  -1 999952824E-01  1 099947900E*ftl 

CHARACTERISTIC  VALUE  DIFFERENCES 

2  927072578E-ft4  -1 610310000E-06  7. 245855850E-O5  -9. 9B40550ME-« 

-1 82875718BE-04  9. 989276O00E-05  -1 858971280E-04  -2  737K3ME-M 

-1 128812368E-03  -1 500e8062BE-B3  -A.  797560W0E-«6  5. 2191538ME-D4 

COMPUTED  COEFFICIENTS 

bk 

.  -1  744821047E-01  9  045283279E-01  3  426584ft37E-01  9  389229542E-W. 

1  9. 375H7004E-W  -3. 402454944E-01  ft  6424ft5451£-ftl  4. 974561162E-01 

-7  652318548E-01  -ft  376485929E-01  -9  386866391E-01  1 39607B933E-81 

COEFFICIENT  DIFFERENCES 

ft  34926S750E-A4  2  79425094BE-04  -6. 382607328E-04  7. 69666802BE-ft4 

2  188920584E-03  -1  774648649E-«3  1  784858868E-03  2  543883813E-03 

-ft  125881640E-04  -5. 147B17»?3E-«?  -1 0859B1871£-«3  2  413A49728E-ft3 


V 
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TABLE  2A 

DEFINING  PARAMETERS.  FUNCTION  VALUES  AND  TRUNCATED 
t>4  =  0.866025404  -  0.51 


CHARACTERISTIC  VALUES  , 

Ak 

4  TPQftftrtfliiflr  fti  4  <wvmw vw.nft  _Q  r-fvvvwvw 

X.  r  JUWWWVoi  X.  O.  k 

.4  rvj  j  MMvwMvww>r,flQ  4  flftflQQQflanr  /u 

1.  &,  JwwSPtt “i  *t.  wwwJC”  X  Wwwwww-TJI 

■j  fvwyvvvyvv-  «u  n  ^wwwww .  rvi  a  fv^wvwwy  Ai 


COEFFICIENTS 

-1  7364817886-01 
9  3969262186-01 
-7. 6604444386-01 


9. 848877530E-81 
-3. 428201438E-81 
-6. 427876100E-81 


FREQUENCIES  m 

LAftiWtfMVMW.flfl 

.  TTWOwrowCTw 

?  'WVWWW- .  /vi 

4'WWWwwMr.ftft 
.  I  VWWWU.IW 

SfWMWWMVWMir.ftft 

1/wwnwvwwwir.fti 


ID  - . 

6+p 

1. 

7  ifttvwwwir.ftfl 

.  lDWWWtTW 

5AAVMWkW  .  flfl 

.  doWwWIt+W 

9-VWWWW  ■  ftn 

JWWWtTW 


FUNCTION  VALUES 


1  420201430E-81 
8  6682540406-01 
-9  3969262106-01 


TRUNCATED 


-3. 876274289E+88 
-2  747461794E+08 
5. 250623797E+00 
2  869146729E+00 
-1  854261256E+00 
-4. 821886543E+88 
1  828485094E+00 
8  262361290E+00 
4. 984842467E+00 
1  369611363E-01 
-2  602917172E+00 
-1 814035377E+80 


2  2482816896+00 
7. 686748895E+80 
2  6842567696+00 
i.  670606168E-81 
-3. 629328817E+00 
-5.  247971842E-01 
5  068968458E+00 
3.  816017954E+80 
-5.  019262483E+00 
-5. 819626286E+80 
-7.  829689 36 1£ -01 
3. 167487943E+00 


J.  OiOUWWWL' w 

-2  740000000E+00 

SAglWMMMWWW  .Aft 

.  tJWWW-Tw 
iWWWWt’w 

4  QeaMAMAPjM 

JL  OJUWWM. 1 W 

J  MMAAteflPtAA 

OiWJWw^TW 

L  ****** .ftft 

8A/-/MWWVW>-  ■  ftrt 

,  <d.0€H*W«*!K3tTOT 
.  AHWUWML  ■  w 

4  Aaaftaoiiflflr  m 
X  wl 

A  A/VWWVMWir  1  Art 
t-  WWUWWL  >W 

.4  {MOfiOfiMfiCAflA 
1  O  lWWwrC  W 


VALUES 


7  (^wwvvvv^  wwt 

5/vwvyvvw^-  wvi 

.  twWWwwC'w 

tiOMMAAArwU 


9.  3969262186-81 

C  MMAMMtf  A1 

3  4202014306-81 


VALUES 

A  AlfMWWWW  .  Art 

C  tWWWOWtMW 

.  wWwwwuTw 

C.  wwvwwv  ■  w 

tr* WVMWMMlf  fti 

UWMWWt T>a 

J/AAMMOArWIQ 

.  OcwwWDW.'W 

5'VWVWWVW  fti 

.  dWWWUUU.*m 

Sfl/flftflMMf.flfl 

.  SwyWENM.’TO 

7  AiftaixmiMif.to 

*>•  wlwvwwt 1  W 

SftiflIWMVWif  iftrt 

.e  nKmMflflTiflfl 

Q.  OAwvwwvim  ■  w 
A4 

.  WWw/w.  vl 
7  i/’fltWKKmr.fWI 
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TABLE  2B 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION 


NUMBER  OF  1TERRTI0NS*  11 


RADIUS  OF  DIFFERENCES  4. 089709056E-88 


COMPUTED  CHARACTERISTIC  VALUES 

-1  ?49999984£-«i 
-L  250600044E-O1 

A  AAAAAAAAir  fti 

”4. 


tAMMMArr.flft 

JQQQQQQQQTTiflfl 

9AOAAAQAAiTiAA 


S008O0828E-02 
-9. 999999700E-02 
-2. 000060O13E-O1 


CHARACTERISTIC  VALUE  DIFFERENCES 
-1  619000000E-09 
4  4348eeeeoE-A9 
-i(02eeeeeeE-« 


jr  f  aaaaaaajm-  aft 

*J.  JjWWWwC 
c.  fOWWWOBI.  w 

Oft 

.  JWWWJWfw 


2  792B08666E-18 
-2  997000000E-A9 
1 3UB00000E-« 


COTWTED  COEFFICIENTS 

-t  73648U33E-81 
9. 396926618E-81 
-7. 660444465E-01 


9. 848077510E-01 
-1  429281343E-01 
427875892E-81 


3. 428281762E-91 
8. 660254154E-81 
-9  396926273E-81 


COEFFICIENT  DIFFERENCES 


-4. 465900000E-06 

TOlwww- 

3ri(W^uww]r  aq 
,  ^JWwwi. 


1  971008000E-69 
-8.  717800800E-09 
•9  nrvYWvr  no 

4.  OfOOWWOtTw 


-I  31858B808E-88 
-1  138800000E-06 
8  317000880E-89 


VALUES 


IMMWEMB 

A  QQQQQQQftigiAO 

110000880(081 


-3  138008008C-89 
5  g^0000Q0Q^ 

-8. 08080O800E-18 


9  396926146E-81 
-4  999999783E-81 
3. 428281503E-81 


6. 441080000E-09 
XCCCOOOO0C  PC 
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TABLE  2C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS3  11 

RADIUS  OF  DIFFERENCES  4. 11221P381E-68 


COHVTED  CHARACTERISTIC  VALUES 


Ak 

-L  7498092UE-01 

1  199993105E+00 

-8. 499958878E-02 

3. 

000094968E+00 

-1  252792343E-01 

3.  999798908E+00 

-9  970799277E-02 

5. 

000177948E+00 

-3. 001535103E-01 

y.  wrowwtTW 

-2  000615370E-01 

1 

099947454E+01 

CHARACTERISTIC  VALUE  DIFFERENCES 

-1  907885100E-05 

6.  894740000E-06 

-4. 112157000E-07 

-9. 

496827000E-05 

2  792343460E-04 

2  010924680E-04 

-2  920072316E-04 

-1 

779475100E-04 

1  535103240E-04 

u.  TnmJjMwCnn 

6. 153698600E-05 

5. 

254SS000E44 

COMPUTED  COEFFICIENTS 

-1  744102864E-01 

9.  821542212E-01 

3.  428315832E-01 

9. 

385152680E-01 

9.  402540961E-01 

-3.  413368353E-01 

8  627195752E-01 

-4. 

969461&49ES1 

-7  662656965E-01 

-6  398963351E-01 

-9  392122933E-01 

3. 

403523720E-01 

COEFFICIENT  DIFFERENCES 

7  621083790E-04 

2  653531781E-03 

-8.  U4401530E-04 

1 

177352994E-03 

-5. 614751220E-04 

-6.  833076510E-04 

3.  305828779E-03 

-L 

053815078E-03 

2  214534740E-04 

-2  897274887E-03 

*4.  803276640E-04 

1. 

667770906E-03 
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TABLE  3A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED  VALUES 

b,  =  -0.866025404  -  0.5i 
4 


CHARACTERISTIC  VALUES 

Ak 

—i  4  'vwwyww .  r*n  _q  p/vwywwir 

-i-  f  vWWwW  ui  i.  “O.  JWWOwvt  vc 

i  ocflflaflflaflf  aj  i  (vvvvvvvvtwm  4  <wwmAAAAf 

A-  &«iwwwvbTi  “■  WWwWiTW  "L  IWWWWL'TjI 

aQQOQOQMr  Ai  Q  rWVVVV^VlAr a  iMMMWMiiwyip  (V. 


3^wwwvvy 
C  (WWVWVWvvr^ 

Li<VWMIA<wm^.<U 


! 


COEFFICIENTS 

-1  736481780E-01 
9. 396926210E-81 
-7. 660444438E-01 

FREQUENCIES  wp 

i  rwvvvwvvv^  wvj 

2/WMWVMVWW  •  ftft 

/WAAwfnw 

3/vwywvwv>r .  >v> 
WVWVUWH.  1  tnJ 

4  7<VMkkMWW  ■ 

.  f  vwvwvw.  *  w 

8/WWWWWVMW>f  ,  ftft 
.  ywwwvm.  ■  w 

Lft/VMWMVWW  ,  /Vi 

UJVWMWLmU 

FUNCTION  VALUES 

-3. 891982761E+80 
-Z  789423297E+80 
5. 128935670E+88 
1  135255271E+88 
-L  866537767E+88 
-4  827465934E+68 
1  OU479539E+80 
8.  2U875056E+68 
4. 771525416E+80 
-3. 328885113E+88 
-1  613126577E+80 
-1  819200816E+80 


9  848077538E-81 
-3.  428291438E-81 
-6.  427876100E-81 

wp+6 

t~VWWWWM-  ■  ftft 

JVWWW^MW 

3  4<HIIMMWWMir.flA 
■  lwwwmnn 

4  4 

•  lwwww^Mn 

5'WWMMWMVM-  .  ftft 
•  tWwArwLH.  I  W 

9vvvwvvwr  ■ 

.  >WUWWM.Mff 

14  'WMVWVW  ■  ft4 
UWUUmilTVI 


7  184134420E+80 
7  083032641E+80 
L  447635953E+88 

—A  ffMffiQftiffiflQ 
KXww>nDtTCTJ 

-3.  068903499E+60 
-1  223838564E-81 
4  876067761E+80 
3. 120731368E+80 

f  ywywu cyyerifta 

T>.  fCTOlDWtW 

1  278369529E+88 
-Z  592876424E-01 
3. 553610963E+08 


3. 420281438E-81 
-8.  668254048E-01 
-9.  396926210E-81 


9. 3969262ME-W. 

ff  aj 

J.  vwwiWWCTU 

3. 420201438E-81 


TRUNCATED  VALUES 


.*>  flOMflfMWVr.ftn 
OvwwwwC'lRJ 

>9  TOflflflflftflflr  i  flft 

t.  f wvwwC»w 

.  U.wvwww’W 

L  138800000E+88 

.4 

X  UWWWWA. 1  W 
Otvvwwwfc'w 

L  010800008E+80 

8PJ/VW WMMflr.flH 

.  CiW«WWAA.Tw 

4.  778808000E+08 
-3.  320000808E+80 

t 

-L  810080000E+80 


iUWWOTW. 1 W 

7- 

Ljjfliwxwmiir.flft 
TtUUUWW.  I W 
4  qqvuuvvvv’^m 

A  WJWWWWCmW 

4. 878000000E+86 

w.  r  tWAAWPPC^w 

1  2788888886+88 

a  e/iflmvwyy  ftj 
L 

A  wWWWWC^W 


V 
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TABLE  3B 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION  VALUES 


NUWER  OF  ITERATION  12 

RADIUS  OF  DIFFERENCES  3.  974336422E-88 


COMPUTED  CHARACTERISTIC  VALUES 


-1  749999976E-01 

4 

A.  CvwwwOt '  W 

-8.  500000306E-82 

7  flflftflftftrtftj  r .  <v> 

-1  249999984E-01 

3. 999999996E+00 

X.  WWWWtLiIi 

e  MWWMWVMir.flft 

v.  WWWwnCTW 

uwwwurXTJl 

Q  flOIVMflflflir,  M 

y.  wwuwvi, ttw 

A  /WMMWVWWW  fti 
“t. 

1  lMXXXWLMn, 

CHARACTERISTIC  VALUE  DIFFERENCES 

-2  352000600E-09 

-s.  1000000066-09 

3. 0571000006-09 

_4  agflflflftflw 

X.  lwwwww._Dy 

-i  5530000006-09 

4. 1400000006-09 

4  4330000006-09 

4  516000000E-09 

r  QdQQQQQQAr  iQ 

u. 

i  /WWMMWWWW  44 

"r.  twwwWJt'll 

J  QMAQQAMr  JA 

JL  wwwwt“lB 

*  * 

COMPUTED  COEFFICIENTS 

! 

-L  736481425E-01 

9. 848077481E-81 

3.  4202016126-01 

9. 396926483E-01 

9  396926191E-01 

-3.  428201077E-01 

-8  66025421 5E -01 

-4  999999072E-01 

-7. 660444628E-81 

-6. 427876192E-01 

-9  3969262106-01 

2  429201421E-01 

COEFFICIENT  DIFFERENCES 

-3. 552500600E-88 

4  9220000066-09 

-18220000006-08 

-2.7253000006-08 

L  8570000006-09 

-3.  5264600666-08 

L  753400000E-08 

-12820000086-06 

L  896206000E-08 

9  227800000E-09 

i  rwvwiwr  11 

>  •  wwwwwi.  XX 

9  1700000606-10 

85 
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TABLE  3C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS2  12 


RADIUS  OF  DIFFERENCES  3.  989296129E-06 


COtfUTED  CHARACTERISTIC  VALUES 

3.  000092925E+B0 
4  999550749E+08 


-1  7485B57B9E-ei 
-1 251056584E-01 
-2.  983978B99E-B1 


1. 19990911 4E +00 
4. 008179343E+0B 
9.  000834954E+08 


-8. 510602756E-02 
-1  003142094E-01 
-1  997906700E-01 


CHARACTERISTIC  VALUE  DIFFERENCES 


-i  494291028E-04 
1 058584240E-04 
-L  682190872E-03 


9  088681000E-05 
-L  793427290E-04 
-8.  349541900E-84 


L  068B74961E-04 
3  142093550E-04 
- Z 093292410E-04 


-9  292490008E-46 
4.  492586900E-04 
2  491780000E-85 


COWUTED  COEFFICIENTS 

-L  738778722E-01 
9  381608229E-01 
-7  636056728E-01 


9  830053900E-B1 
-3.  418689038E-01 
377000067E-81 


3. 438151793E-01 
-8. 652325784E-01 
-9  374728553E-01 


9  393933358E-01 
-4. 997842555E-B1 
3. 421273877E-01 


COEFFICIENT  DIFFERENCES 


2  296941720E-B4 
L  532598893E-03 
-Z  438778223E-03 


1  882362966E-83 
-L  512399768E-04 
-5. 079601381E-83 


-9.  958362660E-04 
-7.  928256368E-94 
-2  2205656 92E -83 


2  9928524 30E -04 
-2 157444718M4 
-1 071S46898E-04 
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TABLE  4A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED 
b4  =  -0.866025404  +  0.5i 


CHflRflCTERISTIC  VRLUES  , 

Ak 

tTr/VWWWlf  Qi  A  VWWV MVW.flfl 

(  JwWWTO.  vl  L  tVWWWU.Tw 

4  'Nrwvwvwy  ftj  4  fVvvwvvvvtr^M 

"i.  vi  WWWUUL’w 

3/wvwww>r  Ai  q  fwwwvwwyir .  /Ml 

.  UWWUKwC’Tl  7  OwwWw)w.Tw 


RrfwvwvMvw  /y-> 
_4  fvwyyvvMy  qa 

“L  pa 

-2.  008080886E-81 


COEFFICIENTS 

-1 736481780E-01 
9  396926210E-01 
-7. 660444430E-01 

FREQUENCIES  wp 

Lmaamamp^q 

2/wwwiimmmwit.m 

3WWWWW  .  /v> 

4~wwvwwy  i  wi 
.  1  .  <  -  ■ 

8/WWWWW  . 

.  yVWWWM.'W 

Limwibimmwiht  .  m 
WTOKmKMU 

FUNCTION  VALUES 

-4. 388366992E+00 
-3  391098418E+00 
4. 187727771E+00 
-1  967826556E+08 
-1  682234702E+00 
-4. 720913679E+00 
5.  827163697E-01 
7. 5635752OLE+00 
3  564087599E+00 
5.  731650938E-01 
-2  450620795E+00 
-1  719677684E+00 


9  84807753flE-0i 
-3  420201430E-01 
-6. 427876100E-01 

V6 

j  mMMMriM 

2.  ftWMMWW  .  ftrt 

.  1  w 

4  a  iwmiihmw*  ■  /V, 

.  llMUMNL^RI 

.  lUMAmju.'C  J 

9~wvwwww ,  /w, 
.  •  w 

4  1  "MWWUir.M 

L  iiwwww*.'!)! 


2  107603510E+00 
7. 024855025E+00 
1  515310314E+00 

-3  851611544E+00 
-3  054850742E+00 
-1. 199074821E-01 
4.  880048139E+00 

3  146831743E+00 
-6. 609237724E+00 

3  277408453E+00 
-2  542912346E-01 
3.  555830719E+08 


3  420201430E-01 
-8. 660254040E-01 
-9. 39692621 0E -81 


TRUNCATED 

A  TtllMBMMr.ftft 

"t.  MmKmC'w 

7  ~vwwvwir  ■  An 

"j.  jjwumoc'w 

4iQMAMMrian 

4  n/MilMWWlf  .M 
X.  AWUHIwCTw 

Dwwwwi.  *; 

_ .  mMMMTWU 

r  dWUUMC'W 

SAAAAAAAQAp^XH 
wwvwwgi.TU 

,  JUMMHM. '  w 

3  560000006E^96 

r«X^^nCa>vl 

iQttUUfiflCiXKI 

-1  710000000E400 


VALUES 


Li  /VV^MWMriAi 


9  396926210E-01 

5  000000000£-^ 
3  420201 438E-01 


VALUES 

a  C4  QOOQQQOE±QGk 

jj00000000g  0^ 

14  iQAAQQQflr  iAA 

3.  2?0000000E*f^ 

3. 
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TABLE  4B 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS  12 
RADIUS  OF  DIFFERENCES  4.  384594668E-88 

COMPUTED  CHARACTERISTIC  VALUES 

Ak 

-1  749999978E-01  1  298000009E+00 

-1 250000817E-A1  3.  99999999SE+88 

3<WMWMWWT_fti  A 

.  wWWwcct^n.  y.  wwwwc^w 

CHARACTERISTIC  VALUE  DIFFERENCES 

-2 991000080E-89  -9. 150000000E-O9  2  927200000E-09  -3  120OOO8OOC-89 

1. 715000000E-A9  5. 490000000E-09  7  114800000E-09  -2  820BMME-49 

2  587000000E-09  1 560000^-09  9  228000080E-10  -2  000000000E-10 

COMPUTED  COEFFICIENTS 

bk 

-1  726481104E-A1  3. 848077497E-01  3. 428281821E-01  9. 2969264UE-U 

9  396926478E-81  -3. 4202W061E-81  -2  660254108E-01  5l  000888267E-01 

-7. 660444617E-01  -6. 427876028E-01  -9. 396926256E-*  1 420281478E-01 

COEFFICIENT  DIFFERENCES 

-6  763880000E-0e  3. 299080008E-09  -3. 913700000E-88  -2  O81780000E-W 

-2  683380000E-A8  -3. 686388000E-88  6. 77W0000BE-W  -2  ttSBMME-tt 

1  871800000E-08  -8. 813800000E-89  4  639008000E-89  -2  9S9BMNE-V 


-8. 500008293E-02  3  000888803E+88 

>1 080080071E-81  5.  eeeeee083E4«e 

A  lUUWUlAAAAf  fti  4  IQAQQAMAP^Ai 

"t.  WWwWX.  vl  i.  IWUUWW.Tpl 
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TABLE  4C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION 


NUMBER  OF  ITERATIONS'  12 


RADIUS  OF  DIFFERENCES  4. 263717385E-06 


COMPUTED  CHARACTERISTIC  VALUES 
Xk 

-L  74B875210E-01  K 
-L  249392160E-O1 
-2  997245095E-01 


L  200383221E+60 
4.  600078669E+00 
9. 000642855E+00 


-8.  496628959E-A2 
-9.  976811942E-82 
-1 999S3€272£-ei 


CHARACTERISTIC  VALUE  DIFFERENCES 
-L 12478995BE-04 
-6  078483500E-05 
-Z  754904950E-84 


-1  8322146ME-04 
-7. 866982000E-O5 
-€.  428553S00E-04 


-2  371041210E-05 
-Z  3988B5798E-4H 
-4.  63728250BE-85 


COMPUTED  COEFFICIENTS 

-1.  727340771E-81 
9  374784685E-01 
-7. 668864915E-01 


9. 831958977E-81 
-3. 414593559E-01 
-6. 399448237E-81 


3  418882556E-01 
-8. 648355125E-01 
-9  382810205E-01 


COEFFICIENT  DIFFERENCES 


-9  141088520E-04 
2.  214152546E-03 
4.  2D4852100E-05 


L  611855279E-83 
-5.  6B7871100E-04 
-Z  842786338E-B3 


L  318873810E-04 
-1  189891550E-03 
-1  4H608468E-63 


VALUES 


Z  999974946E+80 

i  00004  okjehtaHQ 

L  099948277E+81 


8  144070008E-B5 
S.  172342000E-04 


9  390710693E-81 
4. 97548394BE-81 
3  416543817E-91 


6. 215516510E-04 
2  451605900E-83 
1  657612920E-04 
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TABLE  5A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED 

Ag  =  -0.1  +  4.51 


CHARACTERISTIC  VALUES 

Ak 

1.  f  JWwWwC  ill 
tmwOwKC.  tU 

-l  250000000E-01 

COEFFICIENTS  bk 

-L  736481780E-01 
-7.  6604444 30E -01 
9  396926210E-01 

FREQUENCIES  w 

i  *  an 

I  OTOWJWWLtw 

C.  AwWWWXnro 

8<vkw>.v>.vyi(M 

a  CHMXtitXk&CUi M 
1.  C  JWhWwv.mU 

3  0 

OWJ OWWJLmto 

4-wvwwvwv , 

.  ^UWUWW.^ 

FUNCTION  VALUES 

-4.  316014532E+80 
-3  412840693E+08 
-L  687028397E+00 
-4  723241390E+00 
4. 133405575E+00 
6.  967907648E-01 
5.  743445488E-01 
7  537319978E+90 
-2.  454665604E+«i 
-1  721812012E+00 
3.  550468257E+00 
4  129392554E+80 


L'wWkvyyy . 

twWWJWtTw 

4nflrtftA'WinriM> 


9.  848077530E-01 
-6. 427876100E-01 
-3.  420201438E-O1 


4  Tiww^iftonr » rtn 

3.  100000000E+00 

9‘w^yiflflriafl 

1_  lcWWWw.Ma 

4.  100000000E+00 

4-v%^v^v..v>r .  .m 

.  rwwwwxw 


2  270195433E+00 

7.  756090997E+00 
-3. 60U34416E+60 
-5. 128592331E-01 

3.  307169151E+00 

8.  252341996E-01 
5. 09940901 8E +08 
4. 007458598E+80 

-7.  597892780E-01 
3. 178142726E+00 
-3.  583788548E+60 
-4.  401089823E+00 


8CQQQQQQQQC 

_o  r¥VYy¥¥V¥V~  04 
c.  wWXWWTi 

1.  CTWJWwXnJi 


3  420201438E-01 
-9  396926210E-01 
8  668254040E-01 


TRUNCATED 

.  ilTOWWtTW 

-3.  410000000E+00 

>4  /•nflflftWWinfl 

1.  MtoJWW.hw 

4  4  TQQttittCttCiQQ 
.  -WWWX  ul 

S'WVVWVWV*  ftj 

.  f  WRWWwJLnM. 

7  530080000E+00 

d.  <t0wWWWtrW 
1.  r  tlWWTO.TW 

Jfg/WMWlflf  iM 

.  JJwWwt’w 

4  JAftlWWMWtftQ 


INITIAL  VALUE  FOR  UPPER  CLOSE  CHARACTERISTIC  VALUE 

X(6< 1>=  -2  060000008E-02  X<6>  2)=  4.  700000000E+00 


VALUES 


.  WUmKUl.'W 
4  4  flQQQQQQQCjjU 

.L 

4CfWWWVW  i  (Vi 

.  JUWtMXRK.'W 


9.  396926210E-01 
3. 420201430E-01 

Sflftrtflftftftflflf  fli 


VALUES 

Zn^VWVWMlTiflfl 
.  tiUWWW-MW 
7  Tr/wwvw.flfl 
r .  f  JWwwtTw 

Cfwwwwwt  W 

5  4  IWMMWW  .04 
.  Itvwww-ui 

8'VVWWVUW  Qi 

.  Cwwwww  ©x. 

4MAAQQQ QOCxQCt 

.  WeWwWLMW 

7CfvwvMflrtr_fli 

'AWWWW-  wl 

7  i'MnflMflBTiflfl 
J.  If  WIW«WC.TW 
7  Cftftftftflftnnfiflfl 

X  JOWwWwu7!^ 
WWwww. 4  w 
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TABLE  5B 


COMPUTED  DEFINING  PARAMETERS 

FROM  "EXACT"  FUNCTION 

NUMBER  OF  ITERATION  14 

RflOIUS  OF  DIFFERENCES 

2. 

290920 758E -08 

COMPUTED  CHARACTERISTIC  VALUES 

A 

-1  750000008E-01 
-2  999999980E-01 
-1  2500000O6E-01 

Ak 

1  200000005E+00 

9  000000004E+00 

2.  999999999E+08 

0.  jwXJUwJCC  Oc 

-2.  000000011E-01 
-9  999999976E-02 

CHflRflCTERISTIC  VALUE  DIFFERENCES 

8.  020000800E-10  -5.  260000000E-09 

-1.  972000000E-09  -4.  240000000E-09 

6.  200000000E-10  L  160000000E-09 

4  956000006E-10 

4  f¥X¥¥¥¥W¥~  tV) 

-2.  270000000E-10 

COMPUTED  COEFFICIENTS 

bk 

-1  726481404E-01 
-7  660444471E-81 

9  296926221E-01 

K 

9  848077649E-01 
-6.  427875872E-01 
-2.  420201 250E-81 

2.  420201611E-01 
-9  2%92b282E-01 

8.  668252986E-01 

COEFFICIENT  DIFFERENCES 

-2  755600000E-08 

4  125000000E-09 
-1  109600000E-08 

-1  18980U000E-08 
-2.  267400000E-08 
-8.  012000000E-09 

-i.  809500000E-08 

7.  275000000E-09 

J.  WwwITO  oj 

97 


VALUES 


t  an 

-  WWwWWttTw 

A  iQQflQQftQflrwU 

i- 

4.  499999998E+00 


■1  630000000E-09 

2  eeeeeeeeeE-ie 

2  170000000E-09 


9  296926210E-01 
2.  420201452E-01 
5  000000112E-01 


1  500000000E11 

2  188000000E-09 
1  122000000E-08 
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TABLE  5C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION 


NUMBER  OF  ITERATIONS  14 


RADIUS  OF  DIFFERENCES  3.  319S00208E-88 


COMPUTED  CHARACTERISTIC  VALUES 

-1  752347025E-01  * 

-2.  982909368E-81 
-1  252153949E-01 


1  199923297E+00 
9.  00087A285E+O0 
4.  000084482E+00 


-8.  496740820E-82 
-2  002315172E-01 
-9.  969653607E-02 


CHARACTERISTIC  VALUE  DIFFERENCES 


2.  2470258  78E-04 
-1  799963161E-03 
2  152948520E-04 


7.  670301000E-05 
-8.  702845800E-04 
-8. 448212800E-05 


-2  259969820E-05 
2  215171860E-04 
-2.  834629290E-04 


COMPUTED  COEFFICIENTS 


-1.  725825646E-01 
-7  632716826E-81 
9  295062947E-01 


9. 829545296E-01 
~€.  269041342E-01 
-3  421755437E-01 


2.  417477000E-81 
-9  401204511E-01 
8  641210258E-01 


COEFFICIENT  DIFFERENCES 


-6  461227600E-05 
-2.  772840382E-03 
1.  862262720E-04 


8.  532132580E-04 
-5.  883475717E-82 
1.  554007480E-04 


2  724420288E-04 
4.  278200660E-04 
L  904378190E-02 


VALUES 


2  999972497E+00 
1  099951521E+01 
4  499991246E+00 


2.  750266000E-05 

8/TifliftflMr  t%r 

-  WWwwt  w) 


9.  385999003E-01 
2. 402979259E-01 
4  992500852E-01 


1. 092720725E-03 
1.  722217148E-02 
7.  499948500E-04 
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TABLE  6A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED  VALUES 
Ag  =  -0.1  +  4.3i 


CHflRflCTERISTIC  VALUES  \ 
-1  758006000E-01 

i.  KWWWt  ill 
tjWMwW-TJl 


1  200000000E+00 

.  OWWWOWLTW 


_0  CQflQmQQflC.  <y> 
0.  JWwwwCTJc 

QQQQftQQ(MC  ^ 

IQQQQflQQQQC  fli 

WWWWwlTtt 


.  UUMIVUWL'w 

14flQQQOQQQCA(H 
.  IWJWKwWLmJI. 

4WMMMWW  l  /m 

.  jvIMHXKRXK  ’  w 


COEFFICIENTS 

-1  736481780E-01 
-7  660444430E-01 
9  396926210E-01 


9. 848877530E-01 
-6. 427876100E-01 
-3.  428201430E-01 


FREQUENCIES  wp 

1<VMVVliViiW  .  |VA 

wwwmTw 
4  .'WVkWk'Vy  .  flfl 

C-  wwwwt^ 

8<vwvwv^w 

1/wwwvky .  iVj 

DTwwWwCToi 

7  ivwMWiflnr 

S  oWWWWX.tw 

4  wwifln<vir . 

iWJWUlWW.MW 

FUNCTION  VALUES 

-4  326975709E+0U 
-3.  448410174E+00 
-1  692892469E+00 
-4.  726158750E+00 
4  030 32478 7E +00 
3  600966857E+00 
5  622342976E-01 
7  492927000E+80 
-2  459640772E+00 
-1  724485631E+08 
3  948227653E+00 
4  247840727E+00 


wp+6 

L  300000000E+00 
3. 100000000E+00 

9  70QUVWUiriQA 

1.  120000000E+01 

.  IrWiWW-Tw 

4c/wwyjVir  i  Oft 

.  JWWWW-UW 


2  279598703E+00 
7  835005768E+0P 

-3.  593796572E+00 
-5. 095 79901 6E -01 
3. 849863317E+00 
2. 108648817E+00 
5. 1126 304 54E +00 
4  114696285E+00 
-7  536819215E-01 

3  18U06350E+00 
-1.  453948617E+00 
-4  9520595 90E +00 


3  420201430E-01 
-9  396926210E-01 
8  660254040E-01 


9  396926210E-01 
3  420201430E-01 

SflftflflflftflflflT  fti 
WWWWW.  iJl 


TRUNCATED  VALUES 


-4  320080000E+00 
-3  440000000E+00 

tz-fwiflftflftftflriflfl 
O^wwvwu  '  w 

-4  720000000E+00 
4  030000800E+00 

3/vwwvra wv./vi 

WwwwWtM* 

5.  660000000E-01 

f 

ltJW«wW.Tw 

-1  720000000E+00 

.  fta 

.  TTOWWWA.HW 

40iflQQQQQflCiflQ 


0  'vwwwwiriQfl 

7rt«WMUMWWW,flft 

.  Ojvwwvw.tW 

3pwMVM>flftflPiafl 
.  JJWllUWW.Tw 

5/WVWVWMVM*  &l 

.  SW5WWW-TU 

3Ai<vnvkwir.flfl 

UWWWw-Mlro 
i  iiflftnftflflftT.flfl 

7e<vMHMWiftflr  fti 

JWWWWL^l 

7  iftftftrtflftftflf  .flQ 
jgftftftftMflf.ftfl 

i.  lUUWWJW.Tw 
_  A  ftgflflftflflftflT.flQ 
“  JJOwWTOmw 


INITIAL  VALUE  FOP  UPPER  CLOSE  CHARACTERISTIC  VALUE 

X(6.1>=  -2  000000000E-02  x<6.2)=  4  500080008E+00 
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TABLE  6B 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS^  15 
RADIUS  OF  DIFFERENCES  2.  828380666E-08 

COfPUTED  CHARACTERISTIC  VALUES 

-1  750000002E-01  K  L  200000883E+00 

-7  /WWMUVMMW  fti  Q  flflflflftflfiflir  .  flfl 

“i.  J.  WWWWJitTW 

-1  250000019E-01  3.  999999999E+00 

CHARACTERISTIC  VALUE  DIFFERENCES 

1 8000080006-10  -3. 0880000006-09  3  506000000E-10  -1  490000000E-09 

5.  0000000006-12  -2  7 70000000E-09  6  780000000E-10  3. 0000000006-10 

L  917000000E-09  7. 800000000E-10  *1  443000000E-09  1 1500660006-09 

COMPUTED  COEFFICIENTS 

D 

-1. 736481548E-01  k  9. 848077581E-01  3. 4202015926-01  9. 396926202E-01 

-7  660444508E-01  -6. 427875986E-01  -9. 396926259E-01  3. 420201433E-01 

9  396926410E-01  -3. 420201325E-01  8  660253875E-01  5. 0060006646-01 

COEFFICIENT  DIFFERENCES 

-2. 321300000E-08  -5. 053800000E-09  -1 615400000E-08  8. 370000000E-10 

7  845000000E-09  -1  144000000E-08  4  863000000E-09  -2  729000000E-10 

-2  001200000E-08  -1 052500000E-08  1. 652000000E-08  -6. 370000000E-09 


8CQQQflflQ7er  7 

2fwwywv>-»r  Qi  4  4  QQQQQQftflTWU 

.  WWWWffTJl  i.  IUWWwwltJI 

-9.  999999856E-02  4  299999999E+00 


102 


AFWAl-TR-80-31 36 


TABLE  6C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION 


NUMBER  OF  ITERATIONS3  15 

RADIUS  OF  DIFFERENCES  2.  002827512E-08 

COMPUTED  CHARACTERISTIC  VALUES 

X. 

-1  747504865E-01  k  1  200245470E+80 

-2. 987862957E-01  8.  999429818E+00 

-1 250158041E-91  2.  999854478E+00 

CHARACTERISTIC  VALUE  DIFFERENCES 
-2.  495934720E-84 
-1.  213604269E-03 

1.  580405 200E-05 

COMPUTED  COEFFICIENTS 

-1  723584277E-01 
-7  607610209E-01 
9  371556185E-01 

COEFFICIENT  DIFFERENCES 

-2  897502600E-A4 
-5  283422139E-03 

2.  537082472E-03 


-2.  45469E800E-04 
5.  781822B80E-04 
1  455221200E-04 


9  824023985E-01 
-6.  399357E28E-01 
-3.  418899093E-01 


2  4048544E9E-0I 
-2.  85184 7210E-02 
-2  10233S730E-04 


-8.  505172644E-02 
-2. 003951963E-01 
-1  000991056E-01 


5  172643570E-05 
2.  951963360E-04 
9.  910564000E-05 


3.  424241207E-01 
-9.  394357272E-01 
8  661640285E-01 


-4.  039776930E-64 
-2  568938260E-04 
-1.  386245400E-04 


VALUES 


2. 000021437E+00 
1. 100012132E+01 
4  380193887E+00 


•3  142725800E-05 
-1  213162800E-04 
-1.  938869900E-04 


9  390667146E-01 
3.  442706322E-01 
4  975241475E-01 


6  259063920E-04 
-2  258489156E-03 
2  475852470E-02 
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TABLE  7A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED  VALUES 


=  -0. 


CHARACTERISTIC  VALUES 

Xk 

1.  i'JWWWwtiJl  1.  tWWWJTOTW 

_7  fli  a  ±  r*k 

y  OPWWWTOC'TJl  7.  wWwwtTTO 

.4  •gflflfliMfty  a  /w^v>v-y  i  /v* 

1.  cjVCTKWWC^Jl  WWWAW-TW 

COEFFICIENTS  b, 

k 

-1  736481 788E -01  9. 848077530E-01 

-7.  6604444 30E -01  -6  427876100E-01 

9.  396926210E-01  -3  420201430E-01 

FREQUENCIES  wp  wp+g 

.  .VI  4  TflflftVMO/y ,  /** 

WwOwTOtTOT  1.  juwwwOlTW 

-•)  7  > /ww vav^v . 

i.  AWWWWW-’w  w.  1  1 W 

8ft«VkWMVy  .  •-W1  Cl 

.  1'OTwWJWa.  1 7.  f 

•>  iwiaMivw  .  rm  7  iwMvwv-iiy  ■  ai 

4  4  ftVWWVMT .  rw>  A  iiviftflAv^r .  aq 

FUNCTION  VALUES 

-4  332885900E+00  2  284859348E+06 

-3  469381 322E+00  7.  884121558E+00 

-1  695655567E+08  -3.  590400983E+00 

-4  727552407E+00  -5.  880369220E-01 

3  990994621E+08  4  222591S70E+00 

6.  045616763E+00  9  10826081 0E-01 

5  556599425E-01  5. 120071349E+00 

7  466278770E+00  4  183791587E+00 

-2  461998985E+08  -7  508429720E-01 

-1.  725776509E+00  3. 182502995E+00 

6  458849612E+00  4  490988009E+00 

4.  387522930E+00  -5. 4U104025E+06 


+  4.2i 


.0  caQQtMQQQC  7  rtflflflflflflflflr.ftn 

0.  JtwTOwJnit  y  OOWOnmLTm 

o  flftftflIVMMW  04  4  IflftftflflMflT.fti 

C.  wiwwito.  w  i.  i.UUWWWLTUl 

1/VVVVVMVVy  fti  4  'WWMWMW.flft 

.  WWWWCvX  T, 


2  420201430E-01  9  396926210E-01 

-9  396926210E-01  2. 420201430E-01 

errnrzCACAAnc  caa  c  fWWVWWW’  ft-f 

.  OWtJWW.  J.  UUWUWUW.  iH 


TRUNCATED  VALUES 

A  0  'VWMWVMWlTiIVi 

JjWWWW.'W  C.  dOTOWlJlOlfcTW 

_7  i/MflmWMf.  /va  7  /WVMWWWl/Vl 

tWWWWWLTOT  i .  OOuuwuw-Tw3 

— 4  /iwMMflftnflf.flrt  _7  cfwwvwww-infl 

J..  OapwwwwlTw  J.  ^AWWVWL  ’W 

4  Tvywywvwy .  /vi  e  (VVWMVMfty  fli 

t  i  twwww-nKi  0.  wwiwm. Hi 

oft  i  ■vwvMWMWiflfl 

.  931WWTO  t  tw  X6wwww.Tw 

£  fljflftftflftftflf  <ftft  Q  inMflftflflnr  m 

0.  CPrwwut'W-  T  7.  IvwvwUM.  tU 

fli  C  4  WUVIVMW.rtft 

.  UWWWW.TJ1  J.  ltwwww CtW 

7l/(WVWMMflT.flfl  i  J/MIftlWWIftftT./lfl 

.  tOwWAW-TTO  iwwwvwi.  *  w 

‘♦WWWW.tOT  f.  JWWWXwC  ul 

7  a /wvwvwy . /vi 

1.  f  J.wWWWJL’W 

6iCMAWW,flfl  A  JMMMMTiM 

.  ‘fJwWwwK’w  t^WVWWL  1  W 

4vwywwf  inn  je  At ooMMflr.Mi 

.  Xiwwww-’W  J.  “IWWwKtw 


INITIAL  VALUE  FOR  UPPER  aOSE  CHARACTERISTIC  VALUE 

X<6,  i>=  -2  000000086C-02  X<6,2)=  4  380000000E+60 
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TABLE  7B 

COMPUTED  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION  VALUES 


NUWER  OF  ITERATIONS*  14 


RADIUS  OF  DIFFERENCES  4.  418275813E-08 


COMPUTED  CHARACTERISTIC  VALUES 


k 

-L  750000022E-01 

1  200000008E+00 

_0  MOQQMCiC. 

3. 

ftftflflftftrtftV.ftfl 

7  W  A| 

-2.  000000013E-01 

L 

iflflflflflflflflf.AJ 

lWUWUUUL^Vl 

-1  25000001 7E-01 

i  r  .  flft 

-9  999999732E-02 

4. 

199999997E+00 

I 


CHARACTERISTIC  VALUE  DIFERENCES 

1 233000000E-09  -8.  350060000E-09 

z  900000000 e-ii  -s 

L  663800808E-09  -5.  400000000E-10 


S.  425000000E-10 
1.  287000000E-09 

c.  OUWWWWL  U7 


WVWMVW  <y\ 
4L.  twwwwut?? 

3IVWVWVWW-  ift 

.  WWWWW.  10 

3  478000000E-09 


COMPUTED  COEFFICIENTS 

-1.  736481282E-01 
-7.  66044456 7E -01 
*  396926637E-&1 


9. 848077799E-01 
-6.  427875891E-01 
-3.  4202W457E-01 


3  420201680E-01 
-9  396926295E-01 
8.  668253690E-01 


9  396926199E-01 
3. 420201451E-01 
5.  000000281E-01 


COEFFICIENT  DIFFERENCES 


-5  778200000E-08 
1  369906000E-08 
-4  271880000E-08 


-2  686900000E-08 
-2. 088300000E-08 
2.  676000800E-09 


-2.  495200000E-08 
8. 544000000E-09 

7  iMuauaflcjM 


X.  Pet. 

-2.  U4000000E-09 
-2  8L3000000E-08 
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TABLE  7C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION 


NUMBER  OF  ITERATIONS*  14 


RADIUS  OF  DIFFERENCES  4.  425460746E-08 


COMPUTED  CHARACTERISTIC  VALUES 

-1  749991 751E -01 
-2  993761382E-01 
-1. 253883169E-01 


1  200131623E+00 
9. 000250163E+00 
4.  000728998E+00 


-8.  58256592 7E -82 
-2. 00200908iE-0i 
-9.  928138200E-02 


! 


CHARACTERISTIC  VALUE  DIFFERENCES 
-8  249080000E-07 
-6.  238617530E-04 
3  8821 689 TOE -04 


-L  316227406E-04 
-2.  501627800E-04 
-7  389976200E-04 


2.  565927058E-05 
2.  009081120E-04 
-6  186179950E-04 


COMPUTED  COEFFICIENTS 

-1.  740718620E-01 
-7.  637119049E-01 
9  465211470E-01 


9.  838335283E-01 
-6. 413197966E-81 
-3  471542339E-01 


3.  423628101E-01 
-9.  383571265E-01 
8.  601462645E-01 


COEFFICIENT  DIFFERENCES 


4.  2268404 70E-04 
-2.  322538088E-03 
-6  838525966E-02 


9.  742247280E-04 
-1  467812278E-03 
5. 134090889E-03 


-3.  426671420E-04 
-1  335494526E-82 
5.  879139490E-03 
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VALUES 


2.000026430(400 
L  100010780E+01 
4. 19942784 1E+00 


-2. 643022000E-05 
-L  077952000E-04 
5. 721590 300E -04 


9.  385671182E-01 
2. 435831738E-01 
5. 055263387E-01 


1  125502835E-83 
-t  563030795E-03 
-5.  526338720E-03 


V 
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TABLE  8A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED  VALUES 

=  O.I  +  4.1  i 


I 


CHARACTERISTIC  VALUES 

Xk 

1.  I  "  vl  1.  tWWWBAt  *  W 

ivxvwwi.^y  Q  flftfA-Aflflnr  i  rw^ 

-•  UWWWWt'Ul  7. 

ftt  A  i  r¥A 

X  tJwTORWK  wl  *r 


COEFFICIENTS 

-1.  726481780E-81 
-7  660444420E-O1 
9  396926218E-81 

FREQUENCIES  wp 

IflflnAvww  i  .-v% 

waRaWWW.mto 

“i 

d.  l^WUwWW.'fw 

M  WVWliViMW  .  .Ill 

0.  Jw!WwWCT»TW 

.  I'WWWWL’W 

4'WWWWif  iflfl 

cOWwaW-mx1 

FUNCTION  VALUES 

-4  339114189E+O0 
-3  4929942 78E+80 
-1  698316196E+80 
-4  728921079E+0O 
4  022387605E+O0 
8  599607671E+08 
5  486969939E-61 
7  435903825E+80 
-2  464278189E+W 
-1  72782S744E+00 
6  324126963E+80 
2  828841518E+08 


b 


k 


9  848077538E-01 
-6.  427876100E-01 
-3. 420201420E-01 

^p+6 

3WW.wW.tW 

9  7iV^Vii¥¥W  tf¥> 

1  128O00000E+01 

3fWVVWVVW  .  I'trt 

.  AnTAWwCtTW 

4iftfl.WWMflT  .  rtTW 
‘♦WL'VVW^A.  ’  W 


2  290547499E+08 
7. 941878526E+08 
-3.  587173504E+08 
-5.  065550659E-01 
5  887452821E+80 
-6  554092991E+O8 
5  128154321E+00 
4.  267158266E+08 
-7.  481370401E-01 
2  183345658E+08 
5  85565840 7E*00 
-5  046279644E+00 


Q  c^yywww 

2aflQflflflQQflp  vu 
OX 

_4 

vwwwwwv  ill 


3  420201430E-61 
-9.  296926210E-01 
8.  660254040E-01 


.  tA’WyWWL 
4  JMftftflflflflT.Ai 

X-  lwwwmcToi 

4  4  ■wwwvy .  aft 

.  lWWWWM-Tw 


9  396926210E-01 
3. 428281428E-01 

5fWVWWW  fti 

.  UWUWlwtlC"! 


TRUNCATED  VALUES 


-4  330008000E+00 

i  0-^WWWV4>TOT 

-4.  720000000E+O0 

*t.  WWOwWLTtW 
.  J?wwWICtW 

5  400000000E-01 
7  420000000E+00 
-2  460000000E+00 

f  iWWWlTW 

6.ww>fwy>r .  <va 

>.  kWawWItOT 


.  TWWwWXmw 

j. 

C  MQMAAAATJl 4 

5/vwVMwwy.flfl 
4-  CC/WVWWM-  i  Oft, 

v.  JJUUUWUL^ 

5  4  nftftftftflflnr  .Oft 

4v<VMVMW>r.flft 

.  dwwRWW.Tw 

-7  40e000000E-01 

7  4  nfltwww.ftfl 

J.  lOWUWWA.^ 

r  ftHWMwwwinr.ftft 

j. 


INITIAL  VALUE  FOR  UPPER  CLOSE  CHARACTERISTIC  VALUE 
X«6.1>=  -2  000O00000E-02  X(6.2>*  4 


•VWWWW  .  nn 

vWWWWH. '  w 


•^r 
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TABLE  8B 

C0MPUTE0  DEFINING  PARAMETERS  FROM  "EXACT"  FUNCTION 


NUMBER  OF  ITERATIONS*  15 

RADIUS  OF  DIFFERENCES  2.  856678292E-08 

COtfUTED  CHARACTERISTIC  VALUES 
Xk 

-1 7588000UE-61  1  2000000U6E+00 

-3.  000800028E-01  9  000000^0E+00 

-1.  250000154E-01  2.  999999990E+00 

CHARACTERISTIC  VALUE  DIFFERENCES 
1  114000000E-09 
2.  791000000E-09 
L  542600006E-88 

COMPUTED  COEFFICIENTS 

-1.  726481394E-01 
-7  660444544E-01 
9  396926154E-01 

COEFFICIENT  DIFFERENCES 

-3  855200000E-08 
L  141600008E-08 
5  603000000E-09 


-5.  610000008E-09 
2.  00000O000E-11 
1.  036000000E-08 


9.  848077635E-01 
-€  42?8?bt?4E-to. 
-3.  420197904E-01 


-1  551700000E-68 
7.  38300Ct800E-09 
-3. 525939000E-07 


1.  224400000E-09 
-L  450000800E-10 
5.  911000000E-09 


2.  420201591E-01 
-9.  396926213H-&1 
8  660253991E-01 


-1  613800000E-08 
2.  320000000E-10 
4.  910000000E-09 


-8.  508000122E-02 
-1  999999999E-01 
-L  000000059E-01 


VALUES 


7  riftft 

y  wwwWOILmw 


-1  220000000E-09 

SflAAAQAAMT  i  A 

.  wwvwW.  10 

-1  181000000E-08 


9. 296926200E-01 
3  420201397E-01 
4. 999996645E-01 


-9.  028000800E-09 
2  224000000E-09 
2. 25537B000E-07 
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TABLE  8C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION  VALUES 


NUMBER  Of  ITERATIONS^  15 

RADIUS  Cf  DIFFERENCES  2  0202211+56-08 


COMPUTED  CHARACTERISTIC  VALUES 


-1  752805004E-01 
-2  988125622E-01 
-1.  279174132E-01 

Ak 

*  1. 2M160501E+06 

9  001144931E+00 

4.  0620522066+00 

-8.  498616739E-02 
-2.  081353880E-01 
-9.  828075794E-02 

2.  999957827E+08 
1  099998211E+01 
4.  098808774E+00 

CHARACTERISTIC  VALUE  DIFFERENCES 

2.  885002680E-04  -1  6050898806-04 
-1. 186437795E-02  -1  1449211506-02 
-1.  886586848E-03  -2  052385980E-02 

-1  383260930E-05 

1  353880440E-04 
-1  719242859E-03 

4  217287000E-05 
1  789110000E-05 
1  191225970E-03 

COMPUTED  COEFFICIENTS 

b, 

-1  729289544E-01 
-7  627874777E-61 
9.  649118792E-01 

k 

9  841144826E-01 
-6  365813508E-01 
-3.  832885728E-01 

3  420238%3E-01 
-9  383805521E-01 

8.  404534054E-01 

9  389234399E-01 
3  426750860E-01 
5.  413090141E-01 

COEFFICIENT  DIFFERENCES 

-7  272235566E-04 
-2.  256965211E-02 
-2.  521925817E-02 

6  932704270E-04 
-6.  20625 9201E -63 

4  1268429826-02 

-3  753322800E-06 
-1  312068886E-03 

2.  557199863E-02 

7. 691810710E-04 
-6.  549430288E-04 
-4. 130901 405E-02 
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TABLE  9A 

DEFINING  PARAMETERS,  FUNCTION  VALUES  AND  TRUNCATED  VALUES 


CHARACTERISTIC  VALU6S 


-1  7500000806-61 

1.  2000006006+60 

_0  Cflofloa-Ay.  <y> 

7  MttVWWMlTiAi 

WWWWWLhW 

-7  QflQQQQQQQC  CM 
>.  WWWWt  51 

0  » fvi 

i.  WwwWLTJl 

4  iQQQftQQQQCifli 

1.  IWWWwLmJI 

1  iJWJWWt'Ol 

V  lY^VW^iy  t  rtfi 

n.  WWJwlwtw 

.4  ■wwyww  Ai 

I.  WWWW-7J1 

i  QCaQftQflQQPlflQ 

COEFFICIENTS 


-1  7264817886-61 
-7  6684444306-01 
9.  3969262106-01 


9  848877538E-01 
-6  4278761806-01 
-3.  4282014206-01 


FREQUENCIES  wp 

i  0000000006+88 

c. 

3  9080000086+08 

1  0900000066+01 

7  iwmv^r .  .y, 
wwJiWt’w 

4  2000000006+68 

FUNCTION  VALUES 

-4  342356007E+60 
-2  5059422866+80 
-1  6996102416+00 
-4  72959356 3E+00 
4  127276863E+00 

6  741516852E+06 
5.  450584266E-01 

7  419104484E+60 
-2  465389712E+06 
-1  727659556E+90 
7  448711571E+06 

2  6288412896+06 


wp+6 

1  3000000806+00 
3.  1000000886+00 
9  3000008806+00 
1  1200000806+01 

4  4000080006+06 


2.  2935667936+00 
7.  9747182486+00 
-1  5856199896+08 
-5. 058363316E-01 
5. 6646191426+60 
-7  1831856706+06 
5. 132459956E+06 
4  3157582606+00 
-7  48831786 9E -01 
3  1844972406+00 
6  92636828 4E +08 
-4.  8590706486+06 


3  420201439E-01 
-9.  39692621 0E-91 
8.  660254040E-01 


9  396926210E-01 
3  420201430E-01 

5IVMMVVMWW  fti 


TRUNCATED  VALUES 


-A.  7  il  f¥¥¥¥V¥Y~  i  flH 

-  AWWWw.hM 
/•ftflflrfvwvwir.ofl 

-A  TOQQQOQQflCiACI 

4  120000000E+00 
6.  7400000066+00 

c.  dfw¥wwnr  rvj 
3  itwwwWC  vi 

7. 4100006866+00 

.7  ijrf¥¥W¥¥¥~  i  flfl 

-L  7200000006+00 
7.4400000006+00 

7  /“•"uwwwwmr  . 

<L.  ocW.wW.Tw 


7  ftTflftflflMflf.flrt 

r  JiWWOWtHW 

_C  flWMWMMWW  fti 

3.  WWWJW  wl 

5//IWWMW)T.(M\ 

W.'WWWwt  ■  w 

.7  IflfllttMttflfAflfi 

f.  lOtWWI-Tw 

5i  ’MMVjftftftflTiAft 

->  1  vw-ww.  rw! 

-7  4800800006-01 

7  *  'VWWWWVIT .  Art 

iOUWWW-Tw 

O.  'CWWWW.V  W 

-Jt 

^7  WUWWW.’W 


INITIAL  VALUE  FOR  UPPER  CLOSE  CHARACTERISTIC  VALUE 

X'6> 1>-  -2  006000000E-02  X<6,2>=  4  2000000006+00 


AFWAL - TR-80- 31 36 


TABLE  9B 

COMPUTED  DEFINING  P/.-AMETERS  FROM  "EXACT"  FUNCTION  VALUES 


NUMBER  OF  ITERATIONS2  18 

RADIUS  OF  DIFFERENCES  1  118965857E-08 


COMPUTED  CHARACTERISTIC  VALUES 

Xk 

-1 749999999E-01  1 200080600E+00  -8. 500000851E-62  1  O860B888OC+08 

2lvwvwv^cr  Ql  O  QQQQQIWwy.ftft  _ 4  QQQQQQQQ QC_A4  4  iflflttQQQflflCtfli 

-1. 249999581E-01  2. 999999994E+00  -9. 999998924E-02  4. 049999972E+00 

CHARACTERISTIC  VALUE  DIFFERENCES 

-1 450000000E-10  -2  seeeeeeeeE-ie  5.  1260000006-10  -2  oeeeeeeeoc-ie 

4  6000060006-10  5. 4000060006-18  -2  890000000E-10  2  O0O66OQ80C-18 

-4  187000860E-88  6  210000000E-09  -1  076100000E-08  2  814000000E-08 


COMPUTED  COEFFICIENTS 

bk 

-1  736481758E-01  9. 848077525E-01  3  428201472E-01  9. 396926254E-01 

-7  660444437E-01  -6.  427876135E-01  -9  396926283E-01  2. 420201412E-01 

9  396924228E-81  -2  4202141796-61  8  668255942E-81  5. 060012727E-01 

COEFFICIENT  DIFFERENCES 

-2. 240000008E-09  5. 070860000E-10  -4  167000000E-09  -4. 2760006866-09 

6  820000600E-10  3. 4640600006-09  -7  230000000E-10  1  769000000E-09 

L  981960000E-07  1. 274920000E-06  -1  962200000E-07  -1 2726640006-06 


AFWAL-TR-80-31 36 


TABLE  9C 

COMPUTED  DEFINING  PARAMETERS  FROM  TRUNCATED  FUNCTION 


NUMBER  OF  ITERhTHDNS=  16 

RflDIUS  OF  DIFFERENCES  9  623412042E-08 

COMPUTED  CHflRfiCTERISTIC  VALUES 
A . 

-1  750272523E-01  <  1. 200112477E+00 

-2  987588140E-01  9  000772484E+00 

-1  317445485E-01  3.  999254738E +0R 

CHARACTERISTIC  VALUE  DIFFERENCES 
2.  725232900E-O5 
-1  241986824E-03 
6  744548536E-03 

COMPUTED  COEFFICIENTS 

-1.  73582i396E-01 
-7.  626946042E-01 
9  374991959E-01 

COEFFICIENT  DIFFERENCES 

-6. 6028379006-05 
-3.  24982884 7E-03 
2  192425143E-02 


-1. 124774200E-04 
-7.  724841500E-04 
7.  452700300E-04 


9  829918250E-01 
-6.  388I92617E-81 
-1. 078903487E-01 


8. 159280490E-O4 
-4.  748248274E-03 
-2.  34129 7943E-81 


-8.  506619756E-02 
-2  000582121E-01 
-1. 02926104 7E -01 


6.  619755680E-05 
5. 021208900E-05 
2.  926184651E-03 


3  423672743E-01 
-9.  378672626E-01 
8  635735663E-01 


-2.  471212840E-04 
-1.  825358396E-03 
2  451837710E-03 
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V 


VALUES 


2. 000007384E+00 
1  099998201E+01 
4.  056833040E+00 


-7  383640000E-06 
1  798520000E-05 
-6.  823039550E-02 


9  390274240E-01 
3. 424253 723E -01 
2. 669755431E-01 


6  651969940E-04 
-4  052292 726E -04 
2. 338244569E-01 
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Figure  16a.  Graph  of  w(u>),  Ag  =  -0.1  +  4.05i 
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TABLE  9D 

FUNCTION  VALUES  DETERMINED  USING  THE  COMPUTED 
CHARACTERISTIC  VALUES  AND  COEFFICIENTS  FROM  TABLE  9C 


COMPUTED  CHflRflCTERISTIC  VflLUES 


-i.  750272523E-81  k 

1  200112477E+00 

-8.  506619756E-02 

3  000007384E+08 

-2  987580140E-01 

9  000772'484E+08 

-2  000502121E-01 

1  099998201E+01 

-1  317445485E-01 

3.  999254730E+00 

-1  02926104 7E -01 

4. 056833040E+00 

COMPUTED  COEFFICIENTS  k 

K 

-1  735821396E-01 

9  839918250E-01 

3  423672743E-01 

9  390274240E-01 

-7  626946042E-01 

-6. 388393617E-01 

-9  378672626E-01 

3  424253723E-01 

9.  374991959E-01 

-1.  0 7890348 7E-01 

8  635735663E-01 

2  669755431E-01 

FREQUENCIES  u) 

P 

(ji) 

p+6 

A  .-WV-Wi.VW  ■  DJ-i 

1.  OTW*W1  T1CK' 

t.  IwwWWTW 

1  UWT."WVLrW 

0. 

Q  A  ftfl 

7.  -MlwWwlItTOT 

i.. 

A  4 

X.  l£OWWwl“I 

-v  WJWWJIU.TOT 

7  fwvwwwvv ,  M 

cWWwwCtW 

4.  400000000E+00 

FUNCTION  VflLUES  BflSED  ON  COMPUTED  COEFFICIENTS 

AND  COMPUTED  CHflRflCTERISTIC  VflLUES 

TRUE  FUNCTION  VflLUES 

-4  339999982E+00 

2.  289999988E+00 

-4  342356007E+00 

2  293566793E+00 

-3  499999882E+08 

7  969999988E+00 

-3.  505942286E+00 

7  974718348E+00 

-1  689999997E+00 

-3.  579999996E+00 

-1  699610241E+06 

-3. 585619009E+00 

-4  720OO0006E+08 

-5.  000000053E-01 

-4.  729593563E+00 

-5. 068363316E-61 

4  119999996E+00 

5.  660060002E+00 

4. 127276863E+08 

5. 6646191 42E +00 

6  739999997E+08 

-7  179999997E+00 

6.  741516852E+00 

-7. 183185670E+00 

5  400600166E-O1 

5. 129999983E+00 

5.  450584266E-01 

5  132459956E+00 

7  410O00138E+00 

4  309999976E+00 

7.  419104484E+08 

4. 315758268E+00 

-2.  45999999  3E +00 

-7.  399999938E-01 

-2  465389712E+60 

-7  468317969E-01 

_i  TjwwvvMrr  .  /v> 

1.  f  tWWWkl.^OT 

3. 179999994E+00 

-1.  727659556E+08 

3  184497340E+00 

7  440000020E+00 

6  919999722E+00 

7.  448711571E+00 

6  926368284E+00 

2  619999997E+00 

-4.  849999999E+00 

2  628841209E+06 

-4  859070048E+00 
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Figure  17.  Comparison  of  Values  of  w(uj)  for  w(w)  Computed  Us 
Exact  Characteristic  Values  and  Coefficients  and 
w(w)  Computed  Using  the  Characteristic  Values  and 
Coefficients  From  TABLE  9C 
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